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ABSTRACT 

We report on multi-frequency linear polarization monitoring of 15 active 
galactic nuclei containing highly relativistic jets with apparent speeds from ~4 c 
to >40c. The measurements were obtained at optical, 1 mm, and 3 mm wave- 
lengths, and at 7 mm with the Very Long Baseline Array. The data show a wide 
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range in degree of linear polarization among the sources, from <1% to >30%, and 
interday polarization variability in individual sources. The polarization proper- 
ties suggest separation of the sample into three groups with low, intermediate, 
and high variability of polarization in the core at 7 mm : LVP, IVP, and HVP, 
respectively. The groups are partially associated with the common classification 
of active galactic nuclei as radio galaxies and quasars with low optical polariza- 
tion (LVP), BL Lacertae objects (IVP), and highly optically polarized quasars 
(HVP). Our study investigates correlations between total flux, fractional polar- 
ization, and polarization position angle at the different wavelengths. We interpret 
the polarization properties of the sources in the sample through models in which 
weak shocks compress turbulent plasma in the jet. The differences in the orien- 
tation of sources with respect to the observer, jet kinematics, and abundance of 
thermal matter external to the jet near the core can account for the diversity in 
the polarization properties. The results provide strong evidence that the optical 
polarized emission originates in shocks, most likely situated between the 3 mm 
and 7 mm VLBI cores. They also support the idea that the 1 mm core lies at the 
edge of the transition zone between electromagnetically dominated and turbulent 
hydrodynamical sections of the jet. 

Subject headings: galaxies: active — galaxies: quasars: individual (0420-014, 
0528+134, 3C 273, 3C 279, PKS1510-089, 3C 345, CTA102, 3C 454.3) — galaxies: 
BL Lacertae objects: individual (3C 66A, OJ 287, 1803+784, 1823+568, BLLac) 
— galaxies: individual(3C 111, 3C 120) — galaxies: jet — polarization 



Introduction 



Magnetic fields play a prominent role in the physical processes that occur in the jets 
of active galactic nuclei (AGN). The leading model for jet production, acceleration, and 
coUimation involves poloidal magnetic fields that are wound up by the differential rotation of 
a rota.ting disk or ergosphere surrou nding a central supermassive black hole (e.g.. iMcKinney 
20061 : iMeier. Koide. fc Uchidal 120001 ). The twisted field propagates outward as Poynting flux 
in th e polar directions, with eventual conversion into a w ell-focused relativistic plasma flow 



[e.g., IVlahakis fc Konigll |2004J : iMeier fc Nakamura 1120061 ). Within this zone, the magnetic 
field should maintain a tight helical pattern. 

Beyond the jet acceleration region — which may extend over hundreds or tho usands of 
gravitational radii from the black hole (IVlahakis fc Konigll |2004J : iMarscherl |2006| ) — the jet 
may become turbulent or subject to velocity shear. In the former case, the magnetic field 
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should be chaotic, with any hne of sight passing through many turbulent cells and significant 
differences in both strength and direction of the field in adjacer it cells. In co ntrast, velocity 
shear stretches and orders the field lines along the fiow (e.g., iLaing I Il980l ). Shock waves 



passing through the fiow (or vice versa) will compress the component of the field that is 
parallel to the shock front, which imposes order even on a magnetic field that is completely 
chaotic in front of the shock. 

The geometry and degree of order of the magnetic field are therefore key indicators 
of the physical conditions in a jet. Because the primary emission mechanism at radio to 
optical wavelengths is synchrotron radiation, the linear polarization of the continuum can 
be used as a probe of the magnetic field. We can also use the fractional polarization and 
direction of the electric vector position angle (EVPA) to identify distinct features in the jet 
that are observed at different wavelengths. The most prominent features on VLBI images 
of jets in radio- loud AGN are (1) the core, which is the bright, very compact section at the 
narrow end of a one-sided jet, and (2) condensations in the flow that appear as bright knots, 
often called "components" of the jet. The core likely lies some distance from the central 
engine of the AGN, probably either near the end or beyond the zone of acceleration and 
collimation of the jet (jMarscherll2006l ) . The knots usually separate from the core at apparent 
superluminal speeds, but roughly stationary knots are also present in many jets, perhaps 



representing standiri g shocks (e.g.. iJorstad et al.ll200ll : iKellermann et al.ll2004l : iLister 112006 



Jorstad et al. I l2005l . hereafter JOS). The apparent speeds of the components can exceed 
Papp ~ 40c (JOS) that requires the Lorentz factor of the flow exceeds /3app- The leading 
mode l identifles the moving knots as propagating shocks, either t ransverse to the jet axis 
(e.g.. iHughes. AUer. fc Aller I ll98Sl . Il989l : iMarscher &: Gear I llQSSi ) or at an oblique angle 



(IHughes ll200Sl ). Polarization studies can help to determine whether this model is viable. 



The standard paradigm of transverse shocks propagating down a relativistic jet can 
explain many aspects of t he time variability of the brightness , polarization, and structure 



at radio wavelength s (e.g., iHughes. Aller. &: Aller lll98Sl . Il989l : ICawthorne fc Wardle Ill988 



Wardle et al. lll99Sl ). However, this model makes predictions — e.g., that the magnetic fleld 
of a knot should be transverse to jet axis — that often do not match observations. Another 
possibility is that there is a systematicall y ordered component of the rnagnet ic fleld, for 



example one with a helical geometry (e.g., iLyutikov. Pariev. fc Gabuzda ll200Sl ). in the jet 



that modulates its brightness and polarization variability. In fact, a helical magnetic field is 



required in the magnetoh y drodynami c al jet launching models (IMeier. Koide. fc Uchidall2000 



Vlahakis fc Koniglll2004j ). iGabuzda I (120061 ) provides some support for such a geometry in 
recent low frequency VLBI observations of Faraday rotation in AGN jets. However, this 
structure can be related to the magnetic field in the Faraday screen surrounding the jet since 
high frequency radio polarization mapping reveals aspects that are difficult to explain by a 
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helical field flHughes II2OO5I : IZavala fc Tavlor 1120051 ) 



The situation is even more confusing in tlie optical and near-IR regions, where varia- 
tions in flux and pol arization are often extremely rapid, and sub-milliarcsecond resolution 
is not available (e .g., lHagen-Thornlll980l : iMoore. Schmidt, fc Westlll987l : ISmith et al.l 119871 : 
Mead et al.lll990l ). Furthermore, existing optical data ind icate that th e connection between 
variations in brightness and polarization is tenuous (e.g.. ISmith Ill996l ). These factors com- 
plicate interpretation of the optical/near-IR emission, since detailed models are needed but 
the data are not extensive enough to guide them. An alternative approach is to inves- 
tigate the global p olarization behavior across a sizable range of the electromagnetic spec- 



trum. Such studies ( Wills et al. Ill992l : iGabuzda. Sitko. fc Smithlll996l : iLister fc Smith II2OOO 



Gabuzda et al. II2OO6I ) have produced convincing evidence for a connection between the ra- 
dio and optical polarization properties of AGN, suggestive of a common, probably cospatial, 
origin for the emission at these two wavebands. However, such studies have been based on 
single-epoch measurements, whereas the investigation we present here involves multi-epoch 
observations of fifteen objects. This allows us to use both the time and frequency domains 
to explore the geometry and degree of ordering of the magnetic field and other properties of 
the emission regions in relativistic jets. 

Bright jets with high apparent superluminal motion are a prevalent feature of blazars, a 
classification that includes BL Lac objects and optic ally violently variable quasars (OVVs ) 



that are quite rare in the general AGN population (iKellermann et al.l 120041 : ILister II2006I ). 
The prominence of the jets from radio to optical wavelengths and the pronounced variability 
of their polarized emission makes blazars and radio galaxies with blazar-like behavior ideal 
objects for probing the magnetic fields in jets through multi-epoch polarization studies. 

We have carried out a 3-yr monitoring program of fifteen radio-loud, highly variable 
AGN. The program combines roughly bimonthly, high resolution polarized and total inten- 
sity radio images with optical, submillimeter-wave (sub- mm), and millimeter- wave (mm) 
polarization observations performed at many of the same epochs. The sample includes ob- 
jects that are usually brighter than about 2 Jy at 7 mm and 1 Jy at 1 mm, with a mixture of 
quasars (0420-014, 0528+134, 3C 273, 3C 279, 1510-089, 3C 345, CTA 102, and 3C 454.3), 
BL Lac objects (3C 66A, OJ 287, 1803+784, 1823+568, and BL Lac), and radio galaxies 
with bright compact radio jets (3C 111 and 3C 120). In J05, we investigated the properties 
of the jets on milliarcsecond scales and determined the apparent velocities of more than 
100 jet features. For the majority of these components, we derived Doppler factors using a 
new method based on comparison of timescale of decline in flux density with the light-travel 
time across the emitting region. This allowed us to estimate the Lorentz factors as well as 
viewing and opening angles for the jets of all of the sources in our sample. Here we apply 
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these parameters in our interpetation of the multifrequency polarization properties of the 
sources. Because of difficuhies in coordinating the schedules at the different telescopes, our 
observations at different wavelengths are contemporaneous (within two days to two weeks 
of each other) rather than exactly simultaneous. Nevertheless, the program provides the 
richest multi-epoch, multiwaveband polarization dataset compiled to date. 



2. Observations and Data Reduction 

We have measured the total flux densities and polarization of the objects in our sample 
in four spectral regions: at 7 mm (43 GHz), 3 mm (86 GHz), 0.85/1.3 mm (350/230 GHz), 
and optical wavelengths (an effective wavelength of ~600-700 nm). 



2.1. Optical Polarization Observations 



We carried out optical polarization and photometric measurements at several epochs 
using the Two-Holer Polarimeter /Photometer (jSitko. Schmidt, fc Steirull985l ) with the Stew- 
ard Observatory 1.5 m telescope located on Mt. Lemmon, Arizona, and 1.55 m telescope 
on Mt. Bigelow, Arizona. This instrument uses a semi- achromatic half-waveplate spinning 
at 20.65 Hz to modulate incident polarization, and a WoUaston prism to direct orthogonally 
polarized beams to two RCA C31034 GaAs photomultiplier tubes. For all polarization mea- 
surements we used either a 4'.'3 or 8'.'0 circular aperture and, except for 3C 273, no filter. The 
unfiltered observations sample the polarization between ~320 and 890 nm, with an effective 
central wavelength of ~600-700 nm, depending on the spectral shape of the observed object. 
In the case of 3C 273 a Kron-Cousi ns R filter (Xeff ^ 640 nm ) was e mployed to avoid major 



unpolarized emission-line features (jSmith. Schmidt. &: Allen Ill993l ). We accomplished sky 



subtraction by nodding the telescope to a nearby (<30"), blank patch of sky several times 
d uring each observati on. Data reduction of the polarimetry followed the procedure described 



m 



Smith et all (119921 ). 



When conditions were photometric, we obtained differential V^-band photometry {R 
band photometry in the case of 3C 273) of several objects, employing eith er an 8^.^0 or 



16^.^0 circular aperture. We u sed comparison stars in the fields of the AGN ( ISmith et al. 



19851 : ISmith &: Balonek Ill998l ) to calibrate nearly all of the photometry. We determined the 
V magnitudes of 3C 111 and 3C 120 on 1999 Februa ry 13 from the photometric solution 
provided by observations of equatorial standard stars (iLandolt Ill992l ). 



Table [T] lists the results of the optical observations. The columns correspond to (1) 
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source, (2) epoch of observation, (3) filter bandpass of the photometry, (4) apparent magni- 
tude, (5) flux density Jopt, (6) filter bandpass of the polarimetry, where W denotes unfiltered, 
or "white-light" measurements, (7) degree of linear polarization m-opt, and (8) polarization 
position angle Xopt- We have corrected the degree of linear polarizat ion for the statistical 



bias in herent in this positive-definite quantity by using the method of IWardle fc Kronberg 



( ll974j ). Usually this correction is negligible because of the high signal-to-noise (S/N) ratio 



of most measurements. 

Because of the high observed optical polarization levels and high galactic latitudes of 
the majority of the sources, interstellar polarization (ISP) is not a significant concern in most 
cases. Interstellar polarization from dust within the Milky Way Galaxy does appear to be 
a major component of the observed optical polarization for 3C 111 {b = —8.8°) and 3C 120 
(b = —27.4°). A star ~ 50" west of 3C 111 shows very high polarization (m = 3.6% ± 0.2%; 
X = 123° ± 2°); this was used to correct the polarimetry for ISP in the line of sight to the 
radio galaxy. The corrected measurements suggest that the intrinsic polarization of 3C 111 
was typically <3% throughout the monitoring campaign. The average polarization for five 
stars within 8 arcmin of 3C 120 yields m = 1.22% ± 0.06% and x = 98° ± 1°. We have used 
this ISP estimate to correct the observed polarization of 3C 120 for interstellar polarization 
and, as for 3C 111, list the corrected measurements in Table[TJ The corrected values indicate 
that 3C 120 had very low polarization (mopt < 0.5%) throughout the monitoring program. 



Impey. Malkan. &: Tapia I (119891 ) have determined the interstellar polarization along the 
line of sight to 3C 273 to be small, but the low levels of polarization observed for this quasar 
require that this estimate of the ISP be subtracted from the i?-band measurements. Table 
[1] lists the corrected polarization for 3C 273. 



2.2. Submillimeter Polarization Observations 



We performed observations at 1.35 and 0.85 mm with the James Clerk Maxwell Tele- 
scope ( JCMT) locate d on Mauna Kea, Ha\y aii using the Submillir aetre Common User B olome- 
ter Array (SCUBA, [Holland et al. 1119991 ) and its polarimeter (IGreaves et al. II2003I ). The 
initial plan was to observe exclusively at 1.35 mm because the sources are almost always 
brighter, the atmospheric opacity is lower, and the sky is more stable than at 0.85 mm. 
However, failure of the SCUBA filter drum in November 1999 forced a switch to 0.85 mm 
thereafter. The polarizatio n properties of blaza rs tend to be very similar at millimeter and 
submillimeter wavelengths (INartallo et al.lll998l ). so the modest change in wavelength should 
not affect our analysis. For convenience we will refer to the data obtained with the JCMT 
as "1 mm" data. 
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The SCUBA polarimeter consists of a rotating quartz half-waveplate and a fixed ana- 
lyzer. During an observation, the waveplate is stepped through 16 positions, and photometric 
data are taken at each position. One rotation takes ~360 s to complete, and the procedure 
results in a sinusoidally modulated signal from which the Stokes parameters are extracted. 
A typical observation consists of 5-10 complete rotations of the waveplate. We achieved 
fiux calibration in the standard manner with observations of planets or JCMT secondary 
calibrators. We measured the instrumental polarization (~ 1%) during each run by making 
observations of a compact planet, usually Uranus, which is assumed to be unpolarized at 
millimeter /sub millimeter wavelengths. 

We performed the initial (nod compensation, fiat-field, extinction correction and sky 
noise removal, if appropriate) and polarimetric stages of the data reduction using the stan- 
dard SCUBA data reduction packages, SURF and S'/T 0. The Stokes parameters were 
extracted by fitting sinusoids to the data, either half-cycle (8 points) resulting in two esti- 
mates, or full-cycle (16 points) yielding only one estimate but generally giving better results 
with noisy data. We removed spurious measurements by performing a Kolmogorov-Smirnov 
test on the collated data and then calculated the degree and position angle of the polariza- 
tion, corrected for instrumental polarization and parallactic angle. Table [2] lists the data 
with the columns corresponding to (1) source, (2) epoch of observation, (3) wavelength A, 
(4) fiux density /imm, (5) degree of linear polarization mimm, and (6) polarization position 
angle Ximm- As with the optical measurements, the values of mimm have been corrected for 
statistical bias. 



2.3. 3-mm Polarization Observations 

From April 2000 to April 2001 we monitored the sources at 86 GHz using the linear 
polarization system on the the Berkeley-Illinois-Maryland Array (BIMA) in Hat Creek, Cal- 
ifornia during unsubscribed telescope time. The data quality was quite variable, and typical 
total integration times were about 10 minutes per source. We omit from our analysis data 
taken during extremely high atmospheric phase variations or bad opacity conditions. We 
observed planets (Mars, Uranus, and Venus), H II regions (W30H and MWC349), and 3C 84 
as quality checks and for fiux calibration. 

We analyzed the data by mapping all four Stokes parameters for each sideband sepa- 
rately, producing images of total intensity, linearly polarized fiux density, degree of polar- 
ization, and polarization position angle for each source at each epoch. Data having S/N <3 



http://www.jach.hawaii.edu/JCMT/continuum/data/reduce.html 
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for both the total and hnearly polarized intensity images at the central pixel of the source 
were omitted. We determined the total and polarized flux densities by fitting a Gaussian 
brightness distribution to the corresponding image. Table E] contains the results of the ob- 
servations after correction of the degree of polarization for statistical bias: (1) source, (2) 
epoch of observation, (3) fiux density /3mm, (4) degree of linear polarization msmm, and (5) 
polarization position angle Xsmm- 



2.4. 7-mm VLBA Observations 

We performed total and polarized intensity imaging at 43 GHz with the Very Long 
Baseline Array (VLBA) at 17 epochs from 25 March 1998 to 14 April 2001. We describe 
the observations and data reduction in detail in JOS, where the total and polarized intensity 
images are presented. Table H] gives the results of fitting the VLBI core seen in the total and 
polarized intensity images by components with circular Gaussian brightness distributions 
(see JOS). Columns of the table are as follows: (1) source, (2) epoch of observation, (3) 
fiux density in the core /7mm, (4) projected inner jet direction Gjet, (S) degree of hnear 
polarization in the core mrmm, and (6) polarization position angle in the core X7mm- We 
determine the projected inner jet direction from the position of the brightest jet component 
closest to, but at least one synthesized beam width from, the core. 

Table [5] gives the polarization of jet components downstream of the VLBI core that are 
either (i) brighter than or comparable to the core at least at one epoch (/™™p >0.5 /7mm) 
or (ii) have detectable polarization at three or more epochs. These data allow us to probe 
the physics of the strongest disturbances that propagate down the jet. Table O contains: (1) 
source, (2) epoch of observation, (3) designation of the component, which follows that of 
JOS, (4) fiux density in the component Ij'^, (S) distance from the core i?, (6) position angle 
relative to the core 0, (7) degree of linear polarization rrij'^, and polarization position angle 
Xt^, in the component. Values of the degree of polarization in both tables are corrected 
for statistical bias. 



3. Observed Characteristics of Polarization 

Tables [Tl|4] show that the linear polarization for all objects at most wavelengths varies 
significantly both in degree and position angle. Figure [1] represents the relationship between 
the degree of polarization at optical, 1 mm, and 3 mm wavelengths from the whole source 
with respect to the degree of polarization measured in the VLBI core at 7 mm. The data 
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plotted in Figure [T] correspond to the closest observations obtained at opt-7 mm, 1-7 mm, 
and 3-7 mm wavelenertlis for each source. These measurements are simultaneous within 1-3 
day. There is a statistically significant correlation at a signiBcarrce level <=0.05 1 between 
fractional polarization in the VLBI core and degree of polarization at short wavelengths. 
The optical polarization maintains the strongest connection to the polarization in the VLBI 
core (r=0.87, where r is the linear coefficient of correlation). This result confirms the strong 
correlation between the polarization l evel o f the radio core at 7 mm and overall optical 
polarization found by lLister fc Smith I (120001 ). Their sample included quasars with both high 
and low optical polarization. However, it raises the question of why the optical polarization 
shows a better connection to the polarization in the VLBI core than does the polarization 
at 1 mm. 



3.1. Polarization Variability 

We have computed the frac tional polarization variability in dex at each frequency 



using the definition employed by lAUer. AUer. fc Hughes I (l2003al ): 



where m™^^ and m™™ are, respectively, the maximum and minimum fractional polarization 
measured over all epochs at wavelength A, and and (T 771 mill are the corresponding uncer- 

tainties. This approach is especially justified in our case, since all objects were observed over 
the same period of time, although the time intervals are different at different wavelengths 
(~2 yr in the optical, ~1 yr at 3 mm, and ~3 yr at 1 mm and 7 mm). Note that eq. ([1]) 
can produce a negative variability index if the degree of polarization or maximum change in 
polarization do not exceed the uncertainties in the measurements. In this case the data fail 
to reveal polarization variability and we set =0. We introduce the polarization position 
angle variability index, V^: 

T/a _ V rr,\ 

~ 90 ' 



^Throughout the pape r we calculate coefficients of linear correlation and use the method of 
Bowker fc Liebermanlll972l for testing the significance of these coefficients. The hypothesis that there is no 
correlation between two variables, r=0, can be rejected at a significance level e if i = \r j \J\ 



2 degrees of freedom and iV is 



_2, where te_i2;N-i is the percentage point of the t-distribution for iV - 



ie/2;Af- 

the number of observations. 
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where Axa is the observed range of polarization direction and cr^i , cr^2 are the uncertainties 
in the two values of EVPAs that define the range. We treat the 180° ambiguities in EVPAs 
such that Ax\ can not exceed 90°. As in the case of V^, if < then we set = 
and conclude that the observations were unable to measure variability in the polarization 
position angle. 

Figure [2] shows the relationship between the polarization and position angle variability 
indices in the VLBI core and at short wavelengths. A statistically significant correlation 
occurs both between V^p^ and Vf^^ and between V^p^ and V^^^ (r=0.62 and r=0.78, respec- 
tively). This result indicates that changes in the ordering of magnetic fields in the optical 
region and VLBI core are correlated, which suggests that the polarization variability has 
the same origin at the two wavebands. The correlation between the values at 3 mm and in 
the VLBI core might be affected by the region of the jet lying somewhat outside the VLBI 
core that contributes significantly to the 3 mm polarized emission. There is no correlation 
between variability indices at 1 mm and in the VLBI core. 



3.2. A Polarization-Based Classification Scheme 

We use the polarization variability indices to classify the sources in our sample with 
respect to their polarization properties. This classification scheme differs somewhat from 
the classical separation of AGN into radio galaxies, quasars, and BL Lac objects, which we 
have used to discuss kinematics in the parsec-scale jets (JOS). The new categorization reveals 
significant differences in the polarization properties (see §7) among the introduced groups 
that would be diluted in the traditional scheme. 

Figure [3] shows that the fractional polarization and polarization position angle variability 
indices in the 7 mm core are strongly correlated (r=0.83) in a way that produces a clear 
separation of the sources into three groups: LVP - low variability of polarization in the radio 
core, V^jnm <0.2, V^^^^ <0.4; IVP - intermediate variability of polarization, 0.2 < ^7*^111 <0.6, 
0.2 < V^y'^ijn <0.7; and HVP - high variability of polarization in the radio core, Vf^^ >0.6, 
^7mm >0.7. (We note, however, that low electric vector variability indices in the LVP group 
could result from large uncertainties in EVPAs connected with low fractional polarization.) 

The LVP group includes the two radio galaxies 3C 111 and 3C 120 plus the low op- 
tically polarized quasar 3C 273. All three sources display low fractional polarization in 
the optical band and in the 7 mm core. However, in the quasar 3C 273 the optical syn- 
chrotron emission is known to be variable and to have m >10% if dilution from essentially 
unpolarized non-synchrotron components such as the big blue bump is taken into account 
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( llmpey. Malkan. &: Tapia Ill989l ). All three LVP sources are highly polarized and variable at 
1 mm, and the fractional polarization of 3C 273 at 3 mm might be as high as 4%. 

The IVP group consists of four (out of five) BL Lac objects and two highly optically 
polarized quasars, 3C 279 and 3C 345. This group, therefore, includes extremely highly po- 
larized blazars, whose polarization in the optical and at 1 mm can exceed 30% and whose po 



larization at 3 mm and in the VLBI core does not drop below 2-3%. Note that lNartallo et al. 



(119981 ) have found no differences in polarization properties of BL Lac objects and compact 



flat-spectrum quasars at short mm and sub-mm wavelengths. 

The majority of the quasars (five out of eight) and one BL Lac object, OJ 287, form the 
HVP group. The linear polarization of the VLBI core in these sources can be very low — at 
the noise level — or as high as >5%. The polarization variability indices, at 3 mm, 1 mm, 
and optical wavelengths are similar to those for the IVP group, although indices indicate 
more dramatic polarization position angle variability at all wavelengths. We did not observe 
high polarization at short mm wavelengths {msmm >10% and mimm >15%) as in the IVP 
group, although very hig h optical polarizatio n (mopt >20%) might occur in these sources as 



often as in IVP sources (IHagen- Thorn! 1 1 9801 ) . 



Figure m presents the distributions of fractional polarization for all of our measurements, 
separated according to the variability groups. Comparison of the distributions shows that (i) 
the peak of the distribution in the IVP group occurs at a higher percentage polarization than 
for the other two groups, independent of wavelength; and (ii) for a given variability group, 
the peak of the distribution shifts to higher fractional polarization at shorter wavelengths 
except for the LVP sources, where the optical polarization is similar to the polarization in 
the VLBI core. 



4. Total and Polarized Flux Density Spectra 

We have constructed single-epoch total and polarized flux density spectra based on mea- 
surements at different wavelengths obtained within 2 weeks of each other (except 1803+784 
and 1823+568, where, at best, there is nearly a month between observations at different 
wavelengths). We have corrected the optical flux densities shown in Table [1] for galactic 
extinction using values of Ax compiled by the NASA Extragalactic Database. The total flux 
densities at 7 mm are corrected for possible missing flux density in the VLBA images using 
single-dish observations as described in J05. The measurements at 7 mm include the core 
and components within the 1% contours of the peak total intensity at a given epoch (J05), 
/jet = + J™™^. The polarized flux density at 7 mm is integrated over the same VLBA 
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image. Therefore, If,, = ^QU + Uf,„Q,et = Q'^"- + Qr", ^^jet = t/™- + 
mjet = /fet//jet, and Xjet = 0.5 tan-i(f/jet/Qjet), where r"'\ Q™'■^ and t/™"-^ are the Stokes 
parameters of the core, and /^^°™p^ Q™™^, and f/™™P are the Stokes parameters of a given 
polarized jet component. Figure [5] presents the total and polarized spectral energy distri- 
butions (SED), while Table [6] lists the spectral indices (based on total flux densities at 
the three mm wavelengths) and aopt/imm, where 5*1, oc z/""". We calculate the spectral indices 



and using the corresponding polarized flux densities, oc u 



mm "^^"^ "opt/lmm 

Table [6] shows that the total flux density spectra at mm-wavelengths are fiat, indepen- 
dent of the type of source, with the majority of spectral indices having |amm| <0.5. This 
confirms that the mm-wave emission is partially optically thick. The optical-mm spectra of 
the HVP sources are steep, with aopt/imm ^1, while the LVP sources possess much flatter 
optical to 1 mm spectra. In the IVP group aopt/imm depends on whether the source is a 
quasar (aopt/imm ~1) or BL Lac object (aopt/imm ~0.5). The polarized flux density spectra 
at mm wavelengths are fiat or inverted, which indicates an increase of polarized emission with 
frequency at mm wavelengths. Figure [6] reveals a strong correlation between ttopt/imm and 
«opt/imm (r=0.92). However, two sources from the LVP group (3C 120 and 3C 273) deviate 
greatly from the dependence ttopt/imm ~ ctopt/imm (solid line in Fig. [6]), with a significantly 
steeper polarized flux density spectral index than for the total flux density. This suggests 
that at le ast two emission components are present in the optical region, one of which is un- 



polarized. Ilmpey. Malkan. fc Tapia I (119891 ) find an increase in the degree of polarization and 
higher variability of the Stokes parameters in 3C 273 at longer optical wavelengths. Such 
behavior is expected if the optical emission consists of a variable synchrotron component 
plus a relatively static blue unpolarized continuum source, such as the big blue bump. The 
decomposition of the optical synchrotron spectrur n from other optical emission s ources is 



readily seen in the spectropolarimetry of 3C 273 (jSmith. Schmidt, fc Allen Ill993l ). In the 
case of 3C 120, the strong dilution of the optical synchrotron polarization is likely caused by 
host galaxy starlight included within the measurement aperture. 

In the majority of the HVP and IVP sources, the spectral indices ttopt/imm ctopt/imm 
are similar to each other. This suggests that a single synchrotron component is responsible 
for the total flux and polarized continuum from millimeter to optical wavelengths. How- 
ever, three of the HVP sources (0420-014, 1510-089, and 3C 454.3) have a slightly steeper 
polarized flux densit y spectrum relativ e to the total flux density spectrum (this was found 
previously as well by Smith et al. 19881 ). In the IVP sources the opposite trend prevails, ex- 



cept for BL Lac. A steeper polarized spectrum in the quasars indicates the possible presence 
of a non-synchrotron component, such as that observed for 3C 273, although its contribution 
to the total flux should be ^10% relative to the synchrotron emission. A flatter polarized 
spectrum likely rules out any significant contribution by a non-synchrotron component, al- 
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though two (or more) synchrotron components might coexist in the IVP sources. For BL Lac, 
contribution to the optical flux from the host galaxy might help to explain thr slightly steeper 
than the total flux spectral index. 



5. Faraday Rotation in the Inner Jet 



Figure [5] presents measurements of the polarization angle x different frequencies 
obtained within 2 weeks of each other and shows that the direction of mm-wave polarization 
rotates with wavelength. In the HVP and IVP sources the EVPA at 7 mm displayed in the 
figure corresponds to that of the VLBI core only, while in the LVP sources it corresponds 
to Xjet, as described in §4. We attribute this rotation to Faraday rotation by a foregroun d 
screen close to the VLBI core (jZavala fc Taylor I |2004J : lAttridge. Wardle. fc Homan II2005I ). 
We define the rotation measure RM by assuming that at 1 mm we see the unrotated direction 
of the polarization and that the emission from 1 mm to 7 mm propagates through the same 
screen. 



5.1. Faraday Rotation in the HVP and IVP Sources 

Figure [7| shows the result of fitting the dependence of the polarization position angle 
on wavelength by a Faraday rotation law: x = Ximm + RMX^. For 0420—014 our program 
does not contain simultaneous observations at three mm wavelengths necessary to estimate 
RM. We have supplemented our dat a by observations o btained with the VLBA at 7 mm and 



1.3 cm in the program BG073 (see iGomez et al.l 120001 ) . for which 0420-014 and 3C 454.3 



were used as calibrators. For 0420—014 we use our polarization measurement at 1 mm on 
1998 December 13 and EVPAs in the 7 mm and 1.3 cm core on 1998 December 3 from 
BG073. The EVPA at 7 mm is consistent with Ximm at epoch 1998 December 11 from our 
program. For 3C 454.3 and BL Lac we use observations different from those plotted in Fig. 
[5l For 3C 454.3 we use observations at 1 mm and 3 mm from our program (2000 June 6 and 
8, respectively) and 7 mm VLBA observations on 2000 June 9 (BG073). These observations 
are closer to each other than the measurements presented in Figure (2000, November 29 
- December 11) and less affected by the variability of bright component B6 seen on the 
images (J05). We have included the contribution of the component in the polarization at 
7 mm because BQ is brighter and more highly polarized than the core, and located within 
0.2 mas of the core. For BL Lac we use observations carried out over the period 2001 January 
17-28 in order to have simultaneous measurements at three mm wavelengths. 
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Table[7]gives the values of correction for EVPAs at 3 mm and 7 mm along with the intrin- 
sic rotation measure for each object, RM° = RM (1 + ^;)^. Some of the RM values have large 
uncertainties that might be attributed to several problems which our observations possess. 
First, the observations are not completely simultaneous, whereas the p olarization position 



angle can sometimes have short timescales of variability even at 7 mm (ID'Arcangelo et al. 



20071 ). Second, the polarization position angles at 7 mm correspond to those of the VLBI 
core only, while the EVPAs at 1 mm and at 3 mm are from the whole source. Although 
we have tested that, at epochs used for estimates of the rotation measure, jet components 
contribute little to the polarization of the core (except 3C 454.3), their contribution to the 
polarized flux at 1 mm and 3 mm is unknown. Third, the structure of the magnetic field in 
the 7 mm core and in the regions emitting at shorter mm wavelength,- especially at 1 mm,- 
might be different. Therefore, the observed rotation could be an intrinsic property of the 
magnetic field structure. 

To test that the derived rotation measures are reliable, we construct the distributions of 
^X37 = Ixsmm — X7mm| for all pairs of observations simultaneous within 2 weeks for xx before 
and after RM correction (Fig. [H]). Note that the distributions do not include the observations 
used to calculate the rotation measures. This avoids the bias that these observations would 
introduce to the distributions. Ideally, the EVPAs after RM correction should align within 
10°-20°, which corresponds to the uncertainties of our measurements. This is a feature of 
both distributions (uncorrected and corrected) for the IVP sources, reflecting the fact that 
in the IVP sources the uncertanties of EVPAs are comparable to the RM correction values. 
However, EVPAs corrected for RM exhibit a slightly better alignment than AX37 before the 
correction, as expected if the EVPA rotation is caused by a Faraday screen. In the HVP 
sources the RM correction significantly improves alignment between Xsmm and X7mm- 69% of 
uncorrected values of AX37 exceed 20° while only 44% of corrected EVPAs have AX37 > 20°. 
The existence of a few large misalingments between EVPAs a fter the RM correction can be 



attributed to variability of RM in the VLBI core of quasars (jZavala fc Taylor II2004J ). This 
partially affects the analysis performed in §7. 

Table [7] shows that the IVP sources have lower rotation measures, {\RM°\) = (1.1 ± 
0.3) X 10^ rad m'^, than the HVP sources, {\RM°\) = (8.5 ± 6.5) x 10^ rad m^^. Two of 
the HVP blazars, 0528+128 and CTA 102, have EVPA corrections at 7 mm close to 90°. A 
90° flip in EVPA can be caused by the spectral properties o f the source - optic ally thin at 



1 mm but optically thick at 7 mm. However, according to iPacholczyk I (119701 ) the degree 
of polarization should decrease significantly (by a factor of ~ 7) when the optical depth to 
synchrotron self-absorption exceeds unity, while both quasars show only a moderate decrease 
of polarization from 1 mm to 7 mm at the epochs used to calculate RM (mimm = 4.7 ± 1.6% 
and m7mm = 1-8 ± 0.5% for 0528+134; mi^m = 6.5 ± 1.9% and m7mm = 4.0 ± 0.7% for 
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CTA 102). 



In a sample of 40 AGN I Zavala &: Taylor I (120041 ) have determined rotation measures 
using VLBI images at 8-15 GHz that are lower than we derive at mm wavelengths. The 
authors suggest that an external screen "in close proximity to the jet" is the most promising 
candidate for the source of the Faraday rotation. In this context, lower rotation measures 
obtained in the VLBI core at longer wavelengths might reflect a strong decrease in thickness 
of the screen with distance from the central engine. This is expected because the location 
of the core shifts with wavelength owing to optical depth effects. We assume a decreasing 
gradient in electron density of the screen, described by a power law: oc d""^, where d 
is the distance from the central engine. The rotation measure depends on the density and 
magnetic field parallel to the line of sight, integrated along the path through the screen 
to the observer: RM oc f Nf.B\\dl. In the absence of a velocity gradient across the jet, the 
magnetic field along the jet scales as B^, oc and the m agnetic field transverse to the jet 



scales as oc d ^ (IBegelman. Blandford. fc Rees 



1984). If the Faraday screen is at least 



mildly relativistic and has a helical field (iGabuzda II2006I ). B^ provides the main contribution 
to the magnetic field component along the line of sight, therefore, B\\ oc d~^. This leads to 
\RM\ oc d'"" under the approximation that / oc d. The location of the core from t he central 



engine depends on the frequency of the observation u as derived, for example, in iLobanov 
( 19981 ): dcom,u = {B^^F/ 1/)^/'^'', where Bi is the magnetic field at distance 1 pc from the central 
engine, F is a function of the redshift and parameters depending on the jet geometry and 
relativistic electron energy distribution, and the equipartion value /c^ = 1- For a source with 
a fiat mm- wave spectrum dcoT:e,v oc z/~\ which yeilds the dependence of the rotation measure 
observed in the core on the frequency of observation as \RMcore,u\ oc u°-. Therefore, values 
of rotation measure obtained with different sets of frequencies should produce an estimate 
of the parameter a. We assume that the RMo derived from polarization observations of the 
core at frequencies i^i, 1^2, and 1/3, where is the lowest one, yields the intrinsic rotation 
measure at the location (icore,i/3- 

Table [7| contains the published rotation measures obtained in the VLBI core at wave- 
lengths longer than 7 mm, RM^^^^^ and RM^q^^, and derived values of a. For sources 
with RM° available at more than two wavelengths, a is calculated using a least-squares 
method. Table [7] does not show any difference in the values of a between the HVP and 
IVP sources. The average (a) = 1.8 ± 0.5 is similar to the value a = 2 expected for out- 
flow in a spherical or conical wind. This implies that an outflowing sheath wrapped around 
a conically ex panding jet is a reasonab le model for the foreground screen, consistent with 
the finding of I Zavala fc Taylor I (120041 ) . The sheath can result from a mildly relativistic 
outer wind that emanates from the a ccretion disk and confines the inner relativistic jet 
(iGracia. Tsinganos. &: Bogovalovl |2005| ) . According to the RM° values in [TJ the thickness 
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of the sheath is higher in the HVP sources with respect to the IVP sources, a situation 
that should assist in stronger confinement of the jet of HVP blazars. The latter may help 
to explain the finding by JOS that the jet opening angles for quasars are smaller than for 
BL Lac objects. 

We have corrected the values of Xsmm and X7mm in the HVP and IVP sources for Faraday 
rotation by applying the corrections listed in Table [7] to the EVPAs given in Tables |3] and |H 
The adjusted values of the EVPAs are used in the subsequent analysis. 



5.2. Faraday Rotation in the LVP Sources 

In the LVP sources, the core at 7 mm is unpolarized, hence the polarization position 
angle of the inner jet is defined by the polarization of VLBI components within a few mas 
of the core. For the epochs shown in Figure these are components CI and cl at ~4 mas in 
3C 111, components t and u at 1-2 mas in 3C 120, and component B2 at 1 mas in 3C 273. 
Figure [9] shows the dependence of x in 3C 273, which yields a high rotation measure, 
RM = (1.6 ± 0.5) X 10^ rad m~^. The value is consistent wi th the high rotation measure 



obtain ed for component Q6/WQ from the 43/86 GHz images by lAttridge. Wardle. &: Homan 



(120051 ) . RM = 2.1 X 10^ rad m~^. However, a significant decrease in the density of the 
Faraday screen with distance from the core discussed above explains the good alignment 
between the EVPA at 1 mm and Xjet when components are seen farther downstream. This 
occurs for 3C 111 and 3C 120: in 3C 111 ximm = -54° ± 6° and Xjet = -58° ± 10°; in 
3C 120 Ximm = —38° ± 4° and Xjet = —35° ± 6° . The implication is that the magnetic 
field in the 1 mm emission region have the same orientation as that in the jet features a 
few milliarcseconds from the core. Very low polarization in the VLBI core at 7 mm, the 
high rotation measure obtained close to the core in 3C 273, and consistency between the 
direction of polarization in the inner jet and at 1 mm in 3C 111 and 3C 120, all suggest 
that in the LVP sources the 1 mm and 3 mm cores have low polarization as well. The low 
polarization of the core could be the result of either fine-scale turbulence or depolarization 
by a very thick, inhomogeneous f oregro und screen with RM > 5x10^ rad m~^, as suggested 



bv lAttridee. Wardle. fc Homan I (|2005[ ). 



The EVPAs of the LVP sources are not corrected for RM owing to the complexities 
discussed above. We use simplification that, when strong polarized components are detected, 
they have left already the high RM zone. 
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6. Correlation Analysis 

6.1. Comparison of Polarization at Different Wavelengths 

We have calculated linear coefficients of correlation between the polarized flux densi- 
ties in the core at 7 mm and overall polarized flux density measurements at (1) optical, f^j, 
(2) 1 mm, /fy, and (3) 3 mm, /f^. These apply to sources having data from at least three 
essentially simultaneous observations at two wavelengths. (Here we consider the observation 
at two wavelengths "simultaneous" if they are obtained within 2 weeks of each other.) We 
do the same for the correlation coefficients r™ of the the degree of polarization in the core 
at 7 mm and the overall fractional polarization in each of the other wavebands, using the 
same subscript designations. For some sources several measurements have been carried out 
at high frequencies within two weeks of a few VLBA epochs. In these cases, we use only the 
observation that is nearest to the corresponding VLBA epoch. Table M contains coefficients 
of correlation /^y, /f^, and /f^ (columns 2-4) and r^, r^, and (columns 5-7). The number 
of points participating in the computation is indicated in parentheses. The coefficients of 
correlation that are signiflcant at a level e = 0.1 are indicated by bold font. 

The correlation coefficients between optical polarization and polarization in the core 
at 7 mm correspond to neither positive nor inverse correlation despite otherwise simi- 
lar overall polarization characteristics (see §3). Blazars are known to have very short 



timescales of polarization variability at optical wavelengths (<1 day, iHagen-Thornl Il980 



Moore. Schmidt, fc Westlll987l : iMead et al.lll99Cll ). Because of this, we suspect that the lack 



of exact simultaneity and the small number of observations affect the correlation analysis in 
the optical band. 

The statistically signiflcant coefficients of correlation listed in Table [8] indicate a strong 
positive correlation either between the polarized flux in the core at 7 mm and overall source 
at 1 mm and 3 mm, or between the fractional polarization in the core and at shorter mm 
wavelengths. However, the number of sources with such behavior corresponds to 25% and 
45% for the 1-7 mm and 3-7 mm data, respectively. The absence of statistically meaningful 
correlations for other sources might be explained by either small-number statistics, e.g., for 
0420—014 or 1510—089, or by contributions from jet components different from the core 
to the polarized flux at 1 mm and 3 mm. Table [9] gives the coefficients of correlation 
of the sources for which inclusion of the jet components (see §4) changes the signiflcance 
of the correlation. Table [9] shows that including highly polarized components within a 
few milliarcseconds from the core dramatically improves the correlation in 1) all the LVP 
sources, 2) two of the HVP sources (CTA 102 and 3C 454.3), and 3) none of the IVP sources. 
For 3C 279 the correlation coefficient between the 7 mm and 1 mm polarization increases 
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(r^ = 0.51), but does not change its significance substantially. 



6.2. Polarization - Total Flux Connection at Different Wavelengths 

We have computed coefficients of correlation between the degree of polarization and 
total flux density, (listed in Table [TU]) . and between the polarized and total flux densities, 
fx (given in Table [TTl) . The coefficients are calculated for each source that has three or 
more simultaneous measurements of the two quantities at wavelength A (in this case the 
observations are completely simultaneous). Since the JCMT observations were performed at 
1.35 mm and 0.85 mm, we have adjusted flux densities and fractional polarization at 0.85 mm 
to 1.35 mm using the spectral indices a^aia and p rovided in Tabl e P Th e data at 1 mm 



are supplemented by the measurements obtained by iNartallo et al.l (Il998l ) at 1.1 mm for 
sources common to both samples: 0420-014, 0528+134, OJ 287, 3C 273, 3C 279, 3C 345, 
1823+568, BL Lac, and 3C 454.3. Coefficients of correlation r^^^ and /ymm correspond to 
the VLBI core only. 

Table [10] reveals that, in general, the total flux density and percentage polarization 
are not well correlated. However, at 1 mm the one statistically signiflcant coefficient (for 
BL Lac) indicates a decrease of fractional polarization when the source brightens, while at 
optical, 3 mm, and 7 mm wavelengths the few correlation coefficients that are meaningful 
correspond to a positive correlation between the total flux density and degree of polarization. 
An increase of the total flux at mm wavelengths is usually connected with the emergence a 



new superluminal component from the VLBI core (jSavolainen et al.ll2002l ). The component 
might have a different direction of polarization than the core and partially cancel the ob- 
served polarization, thus leading to a poor correlation between the total flux and degree of 
polarization. Table [TTl which lists correlation coefficients between the total and polarized 
flux densities, yields a qualitatively different result. Signiflcant coefficients of correlation 
indicate that the polarized flux density increases as the total flux density rises at all wave- 
lengths. The correlation occurs for 44%, 47%, 46%, and 77% of sources at optical, 1 mm, 
3 mm, and 7 mm wavelengths, respectively. 

Figure [10] shows the dependence between the percentage polarization and total flux 
density at 1 mm and in the core at 7 mm for HVP and IVP sources. For all 12 sources, 
the fractional polarization at maximum flux density in the VLBI core exceeds the average 
degree of polarization (dotted line) while at 1 mm eight sources show the inverse relation: the 
fractional polarization at maximum flux density is lower than the average. The exceptions 
are 3C 279, CTA 102, and 3C 454.3, with all having jet components different from the core 
that appear to contribute to the emission at 1 mm. 1803+784 is another exception - it has 
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the smallest range of flux variability at 1 mm. 

Despite poor correlation between the total flux density and percent polarization, - which 
partially can be affected by complex polarization structure of the mm- wave core region during 
ejection of a new VLBI component - the data support the inference that the fractional 
polarization at 1 mm tends to decrease when a source brightens. In contrast, in the core at 
7 mm the degree of polarization rises along with the total flux density. An increase of both 
the degree of polarization and polarized flux density as the VLBI core brightens implies that 
the flux increase is accompanied by ordering of the magnetic fleld in the emission region. The 
two opposing tendencies seen at 1 mm — lower fractional polarization but higher polarized 
flux as the total flux becomes higher — indicate that the flux increases as emission with weaker 
(but non-zero) polarization becomes prominent. 



7. Polarization Position Angle Behavior 

7.1. Comparison of Polarization Position Angle at Different Wavelengths 

For the IVP and HVP sources we have computed values of deviation between optical 
and mm-wave EVPAs, AxaiA2 = Ixai — XA2I) where Ai is the optical wavelength and A2 is 
1 mm, 3 mm, or 7 mm, for measurements obtained within 2 weeks of each other (recall that 
Ximm corresponds to the EVPA in the VLBI core). Figure [11] presents the distributions of 

Figure [11] shows that the IVP sources possess excellent agreement between polariza- 
tion position angle at different wavelengths: 87% of Axoi (8), 92% of A^os (12), and 90% 
of Axo7 (21) fall into the range 0°-20° (the number in parentheses indicates the number 
of observations). The distribution of Axo? for the IVP sour c es ag rees very well with the 



distribution of |xopt — Xcorejmml obtained by iGabuzda et al. I (12006! ) for 11 BL Lac objects 



and 3C 279. This suggests that in the majority of BL Lac objects (i) alignment between 
the optical polarization position angle and EVPA in the VLBI core at 7 mm is a common 
feature, and (ii) good agreement of the EVPA at different wavelengths from optical to mm 
wavelengths is expected. 

For the HVP sources, the result is qualitatively different: 31% of Axoi (13), 55% of A^os 
(11), and 64% of Axoi (22) are located within 20°. The best agreement is observed between 
optical EVPA and polarization position angle in the core at 7 mm, and between optical 
and 3 mm EVPA, while the distribution of A^oi indicates signiflcantly larger misalignments 
between polarization position angles. 
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In the LVP sources the EVPA at 7 mm corresponds to Xjet, defined by jet components 
within a few milharcseconds of the core. Since the EVPAs are not corrected for rotation 
measure, distributions of AxxiX2 are constructed for all possible pairs of wavelengths (Fig. 
[12]). Figure [T^ shows that for the LVP sources the best agreement in direction of polarization 
occurs between 1 mm and 7 mm as well as between 3 mm and 7 mm. This result, plus a 
good correlation between the fractional polarization light curves at these wavelengths found 
in §6.1, imply that in the LVP sources (i) the mm- wave core has weak polarization and (ii) 
high polarization at the short mm wavelengths (>2-3%) originates in a jet component within 
a few mas of the core. The findings require either (1) negligible Faraday depolarization at 
7 mm and intrinsically low polarization in the core at sub-mm wavelengths, or (2) a highly 
Faraday-thick screen near the core region, causing depolarization of the core at 1 mm. The 
latter is a more suitable explanation due to the high Faraday rotation in the region nearest 
to the core, as discussed in §5.2. 

Comparison of EVPAs obtained at different wavelengths indicates that in the highly 
optically polarized sources (the IVP and HVP blazars) the optical polarization position 
angle tends to align with the EVPA in the VLBI core, while the polarization position angle 
at 1 mm in the HVP sources shows the weakest connection with the optical polarization 
position angle. 

7.2. Direction of Polarization relative to the Jet Axis 

We have compared the polarization position angle, x\j fo^' each polarization measure- 
ment with the position angle of the jet projected on the sky, Gjet (Table H]). 9jet corresponds 
to the nearest VLBA epoch to a given polarization measurement. Figure [T^ gives the distri- 
butions of the offsets between direction of polarization and jet axis, |xa— ©jetl, for polarization 
measured at the different wavelengths. A, and for the different groups of objects. In the case 
of the IVP and HVP sources, X7mm corresponds to the EVPA in the core; in the LVP group 
X7mm is equal to Xjet as defined in §4. The EVPAs at 3 mm and 7 mm for the HVP and IVP 
sources are corrected for RM and the distributions at these wavelengths should be similar 
to that at 1 mm if the RM variability is not significant. It is possible though that the differ- 
ence in the number of observations at the mm wavelengths can introduce some discrepancy 
between the distributions. 

Figure [13] shows obvious alignment of the electric vector with the jet direction in the 
IVP sources: 75% of Xopt (36), 59% of ximm (56), 76% of xsmm (29), and 64% of X7mm (102) 
lie within 20° of the inner jet direction (the numbers in parentheses indicate the number of 
observations at each wavelength), although the widest range of offsets is observed at 1 mm. 
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In contrast, the HVP sources do not show any significant relation between EVPA and jet 
direction, and the values of x\ are distributed almost uniformly from 0° to 90° with respect to 
the jet direction at all wavelengths. An exception is the distribution at 1 mm, which contains 
two peaks at intervals 0°-20° and 80°-90°, suggesting that there is a preference for the electric 
vector to lie along or perpendicular to the jet. The difference between the distributions of the 
IVP and HVP sources is similar to the difference between the distributions of core electric 



vector offsets from the jet direction found by iLister fc Homan I (120051 ) for BL Lac objects and 



quasars, respectively, in the MO JAVE survey at 2 cm (Monitoring of Jets in AGN with VLBA 
Experiments). Although our dichotomy is not based simply on the optical classification of an 
object, we conclude that this difference persists from cm to optical wavelengths and implies 
differences in the magnetic fields and/or processes responsible for the polarized emission in 
these objects. 

In the LVP sources the distributions of \x\~ ©jetl at optical and 1 mm wavelengths do 
not support a connection of the EVPA with the jet direction, although the polarization at 
1 mm seems to avoid being parallel to the jet. The polarization position angle at 3 mm and 
in the inner jet at 7 mm clearly reveals a preferential direction perpendicular to the jet: 89% 
of Xsmm (9) and 62% of X7mm (47) values lie within 70°-90° of the jet axis. Although the 
distribution at 3 mm is dominated by the observations of 3C 273, all three sources make a 
significant contribution in the pronounced peak of the distribution at 7 mm. 

The properties of polarization position angle with respect to the jet direction are distinct 
for each group: good alignment of the electric vector with the jet direction at all wavelengths 
in the IVP sources, chaotic behavior of the electric vector in the HVP blazars independent 
of frequency, and electric vector preferentially transverse to the jet direction in the LVP 
objects at 3 mm and in the inner jet at 7 mm. It is noteworthy that the dichotomy at high 
frequencies is similar to that found for regions farther out in the jet at lower frequencies 



(ICawthorne et al. Ill993l : iLister fc Homan II2005I ). 



8. Interpretation 
8.1. Magnetic Turbulence 



Our data are consistent with the assumption that we observe incoherent synchrotron 
radiation across the wavelength range from optical to mm wavelengths. A uniform magnetic 
field throughout the emission regi on yields the max imum possible degree of polarization. 



^max = (a + l)/(a + 5/3) (e.g., iPacholczyk I Il970l ). which is ~70% at mm wavelength s 
(a ~ 0.5, Table [6]) and ~78% in the optical band (a ~ 1.4; llmpey fc Neugebauer Ill988l ). 
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Such a high degree of polarization has never been observed in a blazar. The polarization can 
be significantly reduced if the emission region is turbulent so that it is composed of many 
cells, each containing a roughly uniform magne tic field that is randomly oriented relative 
to that in other cells (IBurn Ill966l : iBurch Ill979f ) . For N cells within a telescope beam, the 
fractional polarization is given by the equation (IHughes fc Miller Ill99ll ): 



/ViV ± m^ax/V2iV. 



(3) 



The results obtained in §§6-7 indicate that the polarization properties, and therefore 
magnetic field structure, of our sources are comparable within a group but diverse across 
different groups in our classification scheme. We have calculated for each group the mean 
and standard deviation of the maximum degree of polarization over sources measured at 
each wavelength, (m™'^^) ± a™'^^, as well as similar values, (m™"^) ± cr™™, for the minimum 
degree of polarization. The values, given in Table [121 show that both the maximum and 
minimum degree of polarization increase with frequency except for the LVP group, within 
which the optical polarization is similar to the polarization in the VLBI core. At a given 
frequency, the fractional polarization is lowest for the LVP objects and highest for the IVP 
sources, independent of whether the polarization is near the minimum or maximum. We use 
the means and standard deviations given in Table [12] and equation ([3]) to determine whether 
the observed parameters can be derived within a cellular model. Table [12] lists the number 
of cells, N™^{- and N^^, that result in the observed degree of polarization under this model. 

We can calculate the number of turbulent cells for any degree of polarization; how- 
ever, the same number should account for the standard deviation as well. We find that this 
requirement is not satisfied at high-polarization states, where the standard deviation is sig- 
nificantly lower than the model predicts. In low-polarization states, the requirement holds at 
all wavelengths except 1 mm, where the standard deviations for the LVP and HVP sources 
are slightly lower than expected according to the derived number of cells. This suggests that 
the high-polarization states involve an additional process that provides some order to the 
magnetic field, while the low-polarization states can be explained within the cellular model. 

Three interesting consequences follow from analysis of Table [121 First, the observed 
difference in degree of polarization among the groups of objects suggests that the emission 
regions of the LVP and HVP sources contain more turbulent cells than is the case for IVP 
objects - either the cells are finer in the LVP and HVP sources, or the emission region 
is smaller in the IVP blazars. Second, the fact that sources within a group have similar 
fractional polarization at different wavelengths implies that the emitting regions partially 
overlap. Third, the fact that fractional polarization decreases with wavelength suggests that 
the emitting region is larger at longer wavelengths. 
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Faraday rotation by a foreground screen can produce depolarizatio n across a beam that 



also results in frequency dependence of the polarization (IBurn I Il966l ). However, Faraday 



depolarization can not account for the frequency dependence among the groups at optical 
and 1 mm wavelengths since the effect is negligible at such short wavelengths. We assume 
that turbulence and difference in volume of emission is responsible for the frequency depen- 
dence at optical and 1 mm wavelengths but the decrease of polarization at longer mm-waves 
could be caused by Faraday depolarization . To calcula te the frequency dependence of the 
po larization we apply a mo del proposed by iBurn I (Il966l ). with the simplifications suggested 
bv lZavala fc Tavlor I J2004h : 



m{%) = mRM=o|sinc(RMA^)|, 



(4) 



where RM values are rotation measures derived in §5.1. We compare the model with the 
polarization at different wavelengths using the simultaneous data employed for RM esti- 
mates. In this case mpi,M=o corresponds to mimm and is different for each source. Figure 
lH shows the frequency dependence of the polarization according to the observations and 
eq. (jl]). Figure [H] demonstrates that although for some sources (3C 279 and 1803+784) 
the dependence of the polarization on frequency can be caused by Faraday depolarization, 
for the majority of the objects the fractional polarization decreases with wavelength much 
faster than eq. (jl]) predicts. It appears that depolarization from a gradient across the beam 
in a Faraday screen alone is unable to explain the observed frequency dependence of the 
polarization. A difference in the size of the emitting regions and/or in turbulence scale of 
the magnetic field is needed to account for the dependence. 



8.2. Shock Model vs. the Observed Polarization 



A commonly occurring structure that leads to ordering of magnetic fields is a shock, 
which compresses the component of the magnetic field that lies parallel to the shock front. 
The significant dichotomy in degree of polarization at high-polariz ation states amo r ig the 
groups might then be explained by differences in shock strengths. iLister fc Smith I (120001 ) 
proposed a similar explanation to account for differences in polarization properties of low and 
high optically polarized radio-loud quasars (LPRQs and HPQs, respectively) , sugg esting that 
LPRQs have weaker shocks than do HPQs. We use the iHughes fc Miller I (jl99ll ) approach 
to estimate the shock strength for the sources in our sample: 



a + 1 



-r ± i^x — ?7" ; QUI yj , , 

"""^ « + 5/3 2-(l-r/-2)sin2e'' ^ ^ 

where t] = '^ishocked/'^unshocked and n is the number density of the plasma; 9' is the viewing an- 
gle of the jet corrected for relativistic aberration: O' = tan~^[sin9o/(r(cos6o — vT^-T^))], 



-2) sin^ 9' 
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where T is the bulk Lorentz factor and Go is the viewing angle of the jet in the observer's 
frame. The use of the viewing angle with respect to the jet axis in the formula applies to 
a transverse plane-wave shock, and the formula is valid for a shock propagating through a 
plasma with a completely chaotic magnetic field. The maximum possible polarization should 
be seen when the shock is viewed along the plane of compression in the frame of the shock 
(i.e., 9' ~ 90°). Columns 3, 5, 7, and 9 of Table [13] list the values of the estimated shock 
strength t] based on the maximum observed degree of polarization (columns 2, 4, 6, and 8) 
at each wavelength. The last column of the table contains the viewing angle G' in the shock 
frame, as derived from the values of (0o) and (F) obtained for each object in JOS. Because 
we do not have multi-color data in the optical region, the computation of r] adopts aopt/imm 
for the value of Oopt, while at 1 mm, 3 mm, and 7 mm we use the three-point spectral index 
ttmm derived from all three mm wavelengths (see Table [H]). 

The shock strengths listed in Table [13] are similar among the sources and rather low, 
?7 ~ 1.1 — 1.2, and hence cannot be the primary reason behind the diversity in degree of 
polarization. Instead, we can explain this by differences in viewing angle relative to the 
plane of compression, although a very high level of polarization (m ^30%) does require 
much stronger compression, 77 ~ 2 — 3. Note that for 3C 66A JOS find a significant range of 
apparent speeds of superluminal components, from l.Sc to 27c. We use the Lorentz factor and 
viewing angle that corresponds to the average of F and Go derived for the fast components 
only (F=27.8 and 6o=3.3°), because (F) and (0o) averaged over all components (JOS) fails 
to produce polarization as large as that observed for any given value of rj. The latter is 
the case for 1823+S68 as well. However, for 1823+S68 moving components within 1 mas of 
the core were not observed (JOS) and the jet parameters are obtained using kinematics of 
two components at ~2 mas from the core. This implies that F and Go in the VLBI core of 
1823+S68 might be different from those derived in JOS. 

Figure [TJ] shows the dependence between bulk Lorentz factor and polarization variabil- 
ity index. There is a correlation between F and V^p^ (r=0.S3, e = 0.1) and between F and 
^7mm (^"=0.72, e = O.OS). Although V^^^ correlates with F (r=0.42), this correlation is sig- 
nificant only at a level e >0.1S. In the plane-wave shock model, the correlation is expected 
if the minimum degree of polarization corresponds to unshocked plasma and the maximum 
to faster, shocked plasma. In the shock-in-jet model, polarized regions at different wave- 
lengths are partially co-spatial and have essentially the same bulk Lorentz factor. Radiative 
energy losses of the relativistic electrons cause the shorter wavelength emission to occupy 



a smaller fraction of t he full shocked region seen at longer wavelengths (iMarscher fc Gear 



198Sl : iMarscher I [1996]) . Since the shorter wavelength emission is assumed to originate from 
fewer turbulent cells, the degree of polarization and its variability should increase with de- 
creasing wavelength according to equation ([3|). 
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Contrary to what we find at the other wavelengths, the polarization variability index 
at 1 mm does not correlate with the Lorentz factor (r=— 0.12; see Fig. [T5|) . Furthermore, 
according to Table [T21 the number of cells needed in the high-polarization state at 1 mm is 
only slightly larger than, or comparable to, that in the optical region. This conflicts with the 
expectations of the shock model, although iV™^^ at 1 mm is much smaller than at 3 mm and 
7 mm, in keeping with the model. The partial discrepancy suggests that another emission 
component in addition to transverse shocks is prominent at 1 mm. A primary candidate is 
the "true" core, i.e., the bright, narrow end of the jet when observed at a wavelength where 
the emission is completely optically thin (Fig. [T6|) . At wavelengths of 3 mm and longer, 
what appears to be the core is most likely a location at or outside the point where the optical 
depth r(z/) ~ 1. The fact that the optical polarization does show the correlation, as well as 
the best alignment of position angle with the electric vector in the core at 7 mm, implies that 
most of the nonthermal optical emission arises in shocks close to the 7 mm core rather than 
in the "true" core. This suggests that the "true" core does not possess relativistic electrons 
energetic enough to produce the optical synchrotron emission. 



8.3. Magnetic Field Structure in the Jets 



In the shock-in-jet model (e.g., iMarscher Sz Gearlll985l : iHughes. AUer. fc AUer Ill985l ) 



we expect a correlation between total flux and polarization owing to ordering of the magnetic 
field in the shocked region if the quiescent jet has a completely chaotic magnetic field. 
However, if the magnetic field in the unshocked region possesses a component parallel to the 
jet axis, a transverse shock will enhance the total flux but partially cancel the polarization as 
it compresses the turbulent component. In this case the polarized flux will at first decrease 
as the total flux rises. If the shocked emission grows strong enough that its polarized flux 
becomes greater than the initial polarized flux of the ambient jet, the EVPA will flip by 90° 
after a minimum in polarized flux, and the polarized flux will subsequently grow with the 
total flux until the shocked emission weakens. The initial phase of decline in polarized flux 
may easily be missed in observations spaced such as ours. If an oblique shock is responsible 
for an outburst, its effect on the polarization will depend on the angle of the shock front 
relative to the ordered component of the ambient magnetic field, the viewing angle, and 
the bulk Lorentz factor of the jet. It is possible, in some cases, for the net polarization to 
i ncrease, with the EVPA becoming closer to the position angle of the jet axis on the sky 
faughesll2005[ ). 



Figure [TO] shows that the maximum of degree of polarization in the core at 7 mm tends 
to occur when the total flux rises, while at 1 mm the opposite relation is more probable. 
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An increase of the total flux in the VLBI core is usually connected with the emergence of 
a new superluminal component from the optically thick part of the core (see JOS). If these 
components represent transverse shock formation in the jet, they should display electric 
vectors aligned with the jet axis. On the other hand, Ximm at minimum flux might reflect 
the magnetic field direction in the unshocked region. Figure [17] {left panel) displays the 
distributions of misalignment angles between the polarization and inner jet directions at 
times of minimum observed flux at 1 mm. Figure [T71 (middle panel) shows the distributions 
of offsets between X7mm and the local jet direction for each superluminal component listed 
in Table [5] at the point when we can first separate the knot unambiguously from the core 
on the image (i.e., farther downstream than 0.1 mas). Figure [T71 {right panel) presents the 
distribution of degree of polarization in the superluminal knots at their maximum total flux 
when the separation from the core >0.1 mas. Figure [T7] indicates that (1) in the "quiescent 
jet" (prominent at minimum 1 mm flux), the magnetic field aligns with the jet for the LVP 
and HVP sources but is transverse to the jet in the IVP blazars; (2) in the shocked regions 
(superluminal knots) of the LVP objects, the magnetic field is longitudinal or oblique to the 
jet, while for the IVP and HVP sources, its apparent direction is perpendicular to the jet. 

The distributions of offsets between the electric vector and jet direction obtained for the 
IVP sources at different wavelengths (Fig. [T3l) . the distribution of EVPAs in the superluminal 
knots, and the high polarization measured in the knots agree with the cas e of magnetic 



fields lying at oblique angles to the jet axis (ILister. Marscher. &: Gear Ill998l ). (Relativistic 



aberration tends to align the polarization vectors closer to the jet axis, creating the illusion 
that the magnetic fields are closer to transverse than is actually the case in the flow frame 
of the emitting plasma .) This result matches th e findi ng for BL Lac objects reported by 



AUer. AUer. fc Hughes I (j2003bl ) and lAUer et al. I (120051 ) that both the linear polarization 



structure in VLBA images and the multifrequency polarized flux light curves obtained with 
single-dish observations are consistent with the oblique shock hypothesis. Note that the IVP 
sources include most of the BL Lac objects plus those quasars with similar VLBI polarization 
structure. 

In the HVP sources the distribution of offsets between the electric vector and jet di- 
rection at 1 mm is double-peaked (Fig. [T5]) - nearly parallel or perpendicular to the jet. 
According to Figure [TTl at minimum flux at 1 mm the magnetic field has a parallel orien- 
tation, which could be associated with the unshocked plasma. A longitudinal compo nent of 



the m agnetic field in a quiescent jet can arise from a velocity shear across the jet flow (ILaing 



19801 ) . Subtraction of the distribution shown in Figure [T7] from that plotted in Figure [T^ 
produces a distribution for the HVP sources at 1 mm that is quite similar to the distribution 
of the IVP sources. Therefore, the magnetic field in the regions responsible for enhanced 
emission, including superluminal knots, lies at a transverse or oblique angle to the jet axis. 
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in agreement with the shock model. The superluminal components in the HVP sources (Fig. 
[T71) have only a moderate degree of polarization (<10%), indicating that even in the shocked 
region, the magnetic field is fairly turbulent. This is consistent with the high amplitude 
and short timescale of the polarization variability seen at all wavelengths. The existence 
of a longitudinal component of the magnetic field in quiescent states suggests that either 
an increase or a decrease in the degree of polarization is possible when a source brightens, 
depending on the shock strength responsible for the fiux enhancement. This diversity of 
behavior reduces the possibility of observing a correlation between total fiux and fractional 
polarization. 

The general misalignment between the inner jet and the polarization vector in the core at 
7 mm and at optical and 3 mm wavelengths for the HVP blazars (Fig. [T^ fails to reveal any 
connection between the magnetic field and jet structure. In the shock model the distribution 
of electric vector offsets is ex pected to show either a good alignment with the jet direction 



(transverse or oblique shock, iHughes. Aller. fc AUer 



or a bimodal shape (conical shock, ICawthorne fc Cobb 



1985 



Lister. Marscher. fc Gearlll998l ) 



199G 



Cawthorne 1 12006| ) . In the 



conical shock wave model, the on-axis polarization in the shocked region corresponds to the 
upstream magnetic field structure, while the off-axis electric field vectors are oblique to the 
axis and cover a wide range of angles that depend on the shock cone opening angle and the 
viewing angle. A slight asymmetry of the jet cross-sectional structure with respect to the 
axis should lead to relaxation of the dependence of the magnetic field on the jet direction. 
Another factor that can contribute to the spread of x with respect to the jet axis is that, 
despite the high resolution of the 7 mm images, the jet direction at the position where the 
emission occurs might be different from the direction inferred from the images. The latter 
implies that HVP sources possess higher bending of the jets on milliarcsecond scales than do 
IVP sources. Conical shocks might explain the appearance of stationary features in the jet 
( Gomez et al.lll997l : iDaly fc Marscher Ill988l ). which are common in the HVP sources (JOS). 
This model suggests that the core at 7 mm does not always represent the emission from the 
surface where the optical depth r is of order unity. Rather, at any given frequency below the 
turnover in the spectrum, it would be either the r ~ 1 surface or the first standing conical 
shock system in the r < 1 section of the jet, whichever is stronger. If so, the model would 
explain why the separation between the core and c losest superlurn inal knot has a frequency 
dependence in some sources and not in others (see lMarscheiil2006l . and references therein). 



The distributions of position angle offsets between electric vectors and jet direction in 
the LVP sources (Figs. [T31 fTTI) indicate the existence of a persistent longitudinal component 
of the magnetic field in the emission regions at 3 mm, 7 mm, in the quiescent jet, and in 
shocked regions. The modest degree of polarization observed in superluminal components 
of the LVP sources requires a very turbulent magnetic field even in the shocked region. JOS 
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found that LVP sources have the widest jet opening angle in our sample, and the quasar 
3C 273 shows significant variation in the apparent speed of superluminal components across 
the jet. The southern jet components are faster than the northern ones by a factor of 3. This 
suggests that a longitudinal component of the magnetic field might arise from velocity shear 
in the ambient jet, or from complex evolution of the shock emission owing to the formation 
of oblique shocks combi ned with t he eff ects of radiative transfer. The latter is seen in the 



numerical simulations of iHughes I (120051 ). 



Using the parameters of jet components Bl and B2 in the quasar 3C 273 derived in 
JOS and the polarization information given in Table [5], we estimate the distance from the 
core, Az, where a degree of polarization p will be obtained as the result of velocity shear 
in the jet flow. We assume that the magnetic fleld starts out in a turbulent state with cell 
size Ax: Az/AR ^ Tf^AR p/(Ax/?rci)! where Fn is the bulk Lorentz factor at the northern 
side of the jet, AR is the half-width of the jet, and P^ei is the relative speed between the 
northern and southern sides. In the plasma frame on the northern side (subscript "N"), 
Pvei = iPs — /5n)/(1 — /3n/9s), where P = y/1 — T^"^. We assume that the Lorentz factor 
of knot Bl corresponds to the speed at the northern side (Fn = 8.3) and the Lorentz 
factor of knot B2 gives the speed on the southern side (Fs = 13.8). According to Table 
O both components possess a similar maximum polarization (pbi=4.0%, pb2=6.2%) when 
the direction of the electric vector is transverse to the jet axis, so we approximate p ~5%. 
Note that this polarization is achieved in both knots at the same distance from the core, 
~0.8 mas. We use the half-opening and viewing angle of the jet derived in JOS {6 = 1.4° and 
Oo = 6.1°, respectively). This exercise gives a ratio AR/ Ax ~40 and size of a turbulent cell 
Ax ~ 0.001 pc. For the quasar 3C 273 we have obtained several groups of closely spaced 
observations at 3 mm near some VLBA epochs (se e Table [3l). The data reveal a ver y short 



timescale (computed according to the deflnition of iBurbidge. Jones, fc O'Delll Il974l ) of the 
fractional polarization variability at 3 mm, At^ar ~ 1-8 day. This yields the size of the 
region responsible for the variability, Axvar ~ 0.001 pc, which is consistent with the scale of 
turbulence in the magnetic fleld. 

In summary, the sources in our sample show signiflcant variety in the distributions of 
offsets between polarization position angle and jet direction. The observed behavior can be 
explained within the context of shock waves compressing highly turbulent magnetic flelds and 
velocity shear stretching fleld lines along a direction parallel to the axis if the jet structure 
and kinematics on milliarcsecond scales are taken into account. 
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8.4. Possibility of an Electromagnetically Dominated Core at 1 mm 

The polarized emission at 1 mm possesses properties distinguished from the other wave- 
lengths: 1) there is no correlation between fractional polarization variability index and 
Lorentz factor of the jet; 2) the maximum degree of polarization tends to occur when the to- 
tal flux is lower than average; 3) the standard deviation of the mean of the minimum degrees 
of polarization over sources in a variability group is lower than expected in the cellular model 
of magnetic field structure. An inverse relationship between the de gree of polarization and 



total fiux is expected in electromagnetically dominated (ED) jets (ILovelace fc Romanova 



20031 ) if the plasma responsible for a rise in the total fiux has a more chaotic magnetic field 
than that of the ED section. This can be the case even if a shock orders an underlying 
turbulent field in the same sense as does the mostly toroidal field in the ED portion of the 
jet, since the polarization of the ED region should be quite high. 

The upstream ED model of the jet is a promising scenario since it explains the high 
Lorentz factors of the fiow needed to produce the high apparent speeds seen in some blazars. 
Acceleration of the jet to such high values of F is accomplished by conversion of Poynting 
fiux to fiow energy through the Lorentz force that occurs on scales that are much larger 



than the gravitational radius of the central black hole (e.g., IVlahakis II2006I ). The tentative 
anticorrelation between total fiux and degree of polarization at 1 mm and rather positive 
correlation between the values at 3 mm and in the core at 7 mm suggest that the accel- 
eration might end between the VLBl cores at 1 mm and 3 mm. In this case the 1 mm 
emission region should possess a large scale, mostly toroidal, magnetic field component, 
which in the plane of the sky should lie either transverse or longitudinal to the jet direction 



( iLyutikov. Pariev. fc Gabuzda II2005I ). This is consistent with the distributions of EVPAs at 
1 mm with respect to the jet direction in the IVP and HVP sources, in which emission from 
the VLBI core at 1 mm seems to dominate the overall source emission at 1 mm (Fig. [T3|) . 
If the section of the jet seen at 3 mm and 7 mm is beyond the transition point where the 
jet becomes mainly hydrodynamical and turbulent, then the polarization behavior at these 
wavelengths will be similar to that discussed in §§8.1-8.3, with shocks partially ordering the 
randomly-oriented component of the magnetic field. The behavior of the optical polarization 
is more similar to that at 3 mm and 7 mm, and the best agreement between polarization 
position angles for simultaneous observations is observed between the optical EVPA and 
polarization position angle in the core at 7 mm. This suggests that the optical emission 
originates in shocks at or downstream of the transition point, e.g., at or between the VLBI 
cores at 3 mm and 7 mm, and implies that high-energy particle acceleration is more efiicient 
in shocks than in the ED region. 
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9. Summciry 

We have performed multi-frequency, multi-epoch, quasi-simultaneous linear polarization 
observations of a sample of AGN that possess highly relativistic jets and high levels of 
polarization at 1 mm. We separate the sample into three groups according to the pattern of 
variability of the polarization in the VLBI core at 7 mm. This classification appears to be 
connected with physical differences between the objects. However, the sample is small and 
future observations of a larger sample are needed to confirm the findings. The properties of 
the groups are as follows: 

1) The objects with low variabihty of polarization (LVP) in the VLBI core unite two 
radio galaxies with superluminal speeds on parsec scales and a quasar with low observed op- 
tical polarization. The LVP sources have measured optical polarization <2% and radio core 
polarization <1%, and produce strongly polarized radio jet components that are responsible 
for the high polarization observed at 1 mm. The LVP inner jets are most likely subjected 
to very high rotation measures, RM > 5 X 10^ rad m^^ that cause depolarization of the 
core region at mm wavelengths. The optical nonthermal emission is greatly diluted by non- 
synchrotron components, probably associated with the Big Blue Bump. The parsec-scale jets 
are characterized by moderate Lorentz factors, near the low end of the blazar range, and by 
viewing and opening angles that are larger than those of a typical blazar. The distributions 
of EVPAs at 3 mm and in the inner jet at 7 mm with respect to the jet axis reveal a preferred 
longitudinal component of the magnetic field that can be explained by velocity shear. 

2) The sources with intermediate variability of polarization (IVP) in the VLBI core 
include BL Lac objects and quasars with preferred direction of the electric vector at all 
wavelengths close to the direction of the inner jet. They possess very high optical polarization 
(>30%) and relatively highly polarized (>2%) VLBI cores at 7 mm. The polarized and total 
flux density spectra from optical wavelengths to 1 mm suggest that synchrotron components 
with fairly flat spectra, a ^0.5, dominate the total emission at these wavelengths. The 
IVP sources possess a moderate rotation measure, RM < 5 x 10'^ rad m~^, in the core 
that is consistent with the lack of emission lines in the spectra of BL Lac objects owing to a 
relatively low column density of gas. The radio jets of the IVP sources are highly relativistic, 
with r >10. The ahgnment of the polarization with the jet at all wavelengths indicates the 
dominance of a transverse or (in the plasma frame) oblique component of the magnetic field. 

3) The category of objects with highly variable polarization (HVP) in the VLBI core 
consists of the OW quasars and the BL Lac object OJ 287. These display strong variations 
in both degree and polarization position angle. The very high optical polarization seen in 
the IVP sources occurs in the HVP sources, but, in contrast with the IVP blazars, low 
polarization (<1%) at mm wavelengths and in the VLBI core is common. Only a single 
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synchrotron component with a steep spectral index, a ~1.0, is required to produce the 
polarized emission from optical to 1 mm wavelengths, although a non-synchrotron component 
could contribute a small fraction of the total optical emission. The HVP sources possess high 
rotation measures in the VLBI core at 7 mm, RM ~ 3 x 10^ rad m~^, implying that the core 
is surrounded by dense thermal plasma. The radio jets of the HVP sources have Lorentz 
factors at the high-end tail of the blazar distribution. The general misalignment between the 
polarization position angle and jet direction indicates that there is little, if any, connection 
between the two. 

The diversity of polarization properties of the sources in the sample allows us to estab- 
lish connections between parameters at different wavelengths, (i) There is good agreement 
between the optical EVPA and that in the VLBI core at 7 mm for the IVP and HVP blazars. 
(ii) The overall polarization of a source at 1 mm and 3 mm strongly correlates with the de- 
gree of polarization in the inner jet of the LVP sources, (iii) The EVPA at 1 mm and 3 mm 
corresponds to that in the inner jet for the LVP sources, (iv) There is a strong correlation 
between the level of polarization variability in the optical region and in the VLBI core at 
7 mm. (v) The optical and 7 mm polarization variability indices correlate with the bulk 
Lorentz factor of the radio jets, (vi) The degree of polarization decreases with wavelength, 
(vii) The distributions of electric vector offsets with respect to the inner jet direction at op- 
tical and mm wavelengths are similar within each variability group, with distinct differences 
among the three groups. 

In general, the time variability across the various wavelengths follows the expectations 
of models in which rather weak shock waves propagate down a jet containing a turbulent 
magnetic field. After taking into account relativistic aberration, we find that the shock fronts 
generally lie at oblique angles to the jet axis. 

The correlations that we have uncovered tightly link the emission at optical and 3 mm 
wavelengths to the VLBI core at 7 mm. At 1 mm there is a possible anticorrelation between 
total flux density and degree of polarization, but not at 7 mm or 3 mm. This suggests the 
presence of a well-ordered magnetic field in a region that is optically thick at A > 1 mm but 
does not contain electrons with energies sufficiently high to radiate at optical wavelengths. 
During outbursts, the 1 mm polarization is weakened by the somewhat more weakly polarized 
emission from a shock moving through the downstream portion of the jet where the ambient 
magnetic field is chaotic. We suggest that the quasi-steady 1 mm emission arises near the end 
of the section where the jet is accelerated by magnetic forces, a process that requires a tight 
helical field geometry. The acceleration of electrons in this region is apparently inadequate 
to produce optical synchrotron radiation - shocks might be required to do this. Figure [T6] 
presents a sketch of a relativistic jet to illustrate our interpretation. 
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Our study demonstrates the utility of multiwaveband polarization monitoring of blazars 
when combined with multi-epoch imaging of their jets with VLBI at mm wavelengths. Fur- 
ther progress is possible through more intense time and frequency coverage of a large number 
of objects. 
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Table 1. Optical Polarization Data 



Source 


Epoch 


Filts 


Mag 


/opt [mJy] 


Filt^ 


mopt[%] 


Xopt[°] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


3C 66A 


1999/02/12 








W 


13.8±0.2 


37.2±0.3 




1999/02/13 


V 


14.70±0.02 


4.81±0.08 


w 


13.3±0.2 


45.2±0.4 




1999/02/14 


V 


14.74±0.02 


4.63±0.08 


w 


14.5±0.2 


47.4±0.3 




2000/09/27 


V 


14.78±G.02 


4.46±0.G8 


w 


29.7±G.2 


17.2±G.2 




2000/09/28 


V 


14.81±0.02 


4.35±0.09 


w 


26.8±0.3 


17.1±G.4 




2000/09/29 








w 


24.2±0.3 


13.4±0.4 




2000/11/30 


V 


14.40±0.02 


6.31±0.11 


w 


23.7±G.l 


25.4±G.l 




2000/12/01 








w 


22.8±G.2 


25.1±G.2 




2001/01/19 


V 


14.43±0.02 


6.16±0.10 


w 


24.6±G.l 


IG.liG.l 


3C 111 


1999/02/13 


V 


18.95±0.12 


0.10±0.02 


w 


3.1±0.8 


126±4 




2000/04/03 








w 


G.8±G.8 


1G6±5 




2000/09/27 








w 


1.5±G.6 


1G4±3 




2000/11/30 








w 


1.9±G.5 


131±3 




2000/12/01 








w 


1.6±0.8 


110±5 




2001/01/19 








w 


1.3±0.4 


89±3 



Note. — 



The table appears in full form electronically. 
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Tablc 2. JCMT Polarization Data 



Source 


Epoch 


A [mm] 


/lmm[Jy] 


TOlmm[%] 


Xlmm[°] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


3C 66A 


1998/05/15 


1.35 


0.38±0.04 


36±4 


24±2 




1998/07/17 


1.35 


0.41±0.04 


19±4 


176±6 




1998/09/30 


1.35 


0.45±0.04 


14±2 


56±4 






J. .ou 


f) 41 _(_n nA 
11.^1 1 111 1^ 




("I 79-1-1 fi'i 




1998/12/13 


1.35 


0.47±0.05 


14±2 


38±4 




1999/02/18 


1.35 


0.43±0.04 


11±2 


157±4 




1999/04/22 


1.35 


0.48±0.04 


9.3±2.5 


134±7 




1999/09/23 


1.35 


0.59±0.06 


9.7±2.5 


158±10 




2000/06/08 


0.85 


0.61±0.06 


9.6±1.3 


9±4 




2001/01/17 


0.85 


0.61±0.06 


11±2 


23±4 




2001/01/23 


0.85 


0.59±0.06 


14±3 


4±5 


3C 111 


1998/07/17 


1.35 


1.03±0.10 


12±2 


39±4 




1998/09/30 


1.35 


1.45±0.10 


3.9±0.7 


117±5 




1998/11/23 


1.35 


1.26±0.13 


4.2±0.9 


123±6 




1998/12/13 


1.35 


1.37±0.14 


5.8±1.6 


39±7 




1999/02/18 


1.35 


1.04±0.10 


4.5±1.1 


126±6 




1999/09/23 


1.35 


0.89±0.09 


6.4±2.0 


143±11 




1999/12/24 


0.85 


0.93±0.10 


(4.9±1.9) 


(17±10) 




2000/06/08 


0.85 


1.05±0.11 


(1.0±1.5) 


(144±39) 




2001/01/17 


0.85 


3.61±0.36 


2.2±0.7 


46±8 




2()()l/()l/23 


0.85 


r).99±().(i0 


i.oio.;-! 


17±7 



Note. — The table appears in full form electronically. 
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Table 3. BIMA Polarization Data 



Source 

(1) 


Epoch 

(2) 


hmm [ Jy] 


(4) 


XSmni [ ] 

(•"') 


3C 66A 


2000/04/21 


0.82±0.01 


3.5± 0.7 


171± 5 




2000/04/25 


0.98±0.01 


5.0± 0.5 


1± 2 




2001/02/06 


0.89±0.01 


9.3± 0.3 


17± 1 




2001/03/20 


0.65±0.01 


4.2± 0.2 


2± 1 




2001/03/29 


0.85±0.01 


5.9± 0.6 


10± 2 




2001/04/03 


0.86±0.01 


5.0± 0.5 


8± 2 




2001/04/07 


0.71±0.01 


4.8± 0.4 


17± 2 




2001/04/12 


0.90±0.01 


3.3± 0.8 


32± 6 




2001/04/16 


0.71±0.01 


3.1± 0.6 


0± 5 




2001/04/20 


1.00±0.02 


5.6± 0.8 


15± 4 


3C 111 


2000/04/04 


2.63±0.01 


0.3± 0.1 


134± 9 




2000/04/14 


3.33±0.02 


2.0± 1.2 


45±14 




2000/12/06 


4.10±0.01 


0.5± 0.2 


138±11 



Note. — 



The table appears in full form electronically. 
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Table 4. VLB A Polarization Data for the Core 



Source 


Epoch 


-^7111111 [ Jy] 




TO7mm[%] 


X7mm[°] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


3C 66A 


1998/03/25 


0.52±0.05 


-146±1 


2.3±0.6 


26±7 




1998/05/30 


0.45±0.07 


-155±1 


3.7±0.8 


13±5 




1998/07/31 


0.30±0.03 


-150.7±0.5 


2.9±0.6 


35±6 




1998/10/05 


0.31±0.04 


-155±1 


2.7±1.1 


36±10 




1998/12/10 


0.39±0.12 


-159±2 


3.1±1.3 


13±11 




1999/02/11 


0.48±0.05 


-146±2 


2.2±0.5 


16±6 




1999/04/29 


0.50±0.03 


-167±1 


1.4±0.6 


-6±11 




1999/07/18 


0.47±0.02 


-162±2 


2.3±1.0 


27±11 




1999/10/06 


0.69±0.03 


-154.4±0.5 


3.2±0.4 


18±4 




1999/12/05 


0.75±0.03 


-159±1 


7.7±0.4 


5±3 




2000/01/24 


0.61±0.05 


-155.4±0.5 


3.9±0.4 


0±3 




2000/04/05 


0.60±0.05 


-158±1 


3.8±0.5 


14±4 




2000/07/17 


0.63±0.08 


-163±1 


4.3±1.7 


28±11 




2000/10/01 


0.73±0.04 


-159±1 


6.1±0.3 


28±3 




2000/12/11 


0.60±0.06 


-160±2 


5.7±1.0 


27±10 




2001/01/28 


0.63±0.07 


-162±2 


6.1±1.5 


35±8 




2001/04/14 


0.60±0.05 


-160±1 


5.4±0.3 


28±3 



Note. — 



The table appears in full form electronically. 
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Table 5. VLBA Polarization Data for the Jet Component 



Source 


Epoch 


Comp. 




i?[nias] 


e[1 




comp ro] 
A. 7mm L J 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


3C 111 


1998/03/25 


CI 


0.43±0.08 


2.50±0.10 


66.8±0.5 


6.3±2.1 


-44±10 






cl 


0.14±0.06 


2.15±0.08 


67.6±0.5 


15.7±2.8 


-43±8 




1998/05/30 


CI 


0.54±0.08 


2.71±0.10 


67±1 


6.8±1.8 


-67±7 






cl 


0.10±0.04 


2.21±0.10 


69±1 


36±6 


-65±6 




1998/07/31 


Cl 


0.32±0.08 


2.93±0.10 


67.7±0.5 


9.4±2.1 


-55±7 






cl 


0.16±0.05 


2.34±0.15 


68.5±0.5 


13±3 


-47±8 




1998/10/05 


Cl 


0.22±0.06 


3.31±0.15 


67.5±1.0 


9.3±2.5 


-45±8 






cl 


0.09±0.04 


2.64±0.15 


70±2 


32±6 


-45±6 




1998/12/10 


Cl 


0.17±0.05 


3.73±0.15 


67±1 


15±4 


-40±10 






cl 


0.16±0.04 


3.16±0.15 


69±2 


8.8±2.1 


-54±10 




1999/02/11 


Cl 


0.17±0.05 


4.01±0.15 


67±1 


18.8±4.5 


-57±10 






cl 


0.15±0.05 


3.41±0.15 


69±2 


8.8±2.1 


-60±10 




1999/04/29 


Cl 


0.12±0.05 


4.33±0.20 


66±2 


21±5 


-22±20 






cl 


0.16±0.07 


3.62±0.20 


70±2 


15.6±4.3 


-41±20 




1999/07/18 


Cl 


0.20±0.05 


4.69±0.20 


66±2 


6.0±1.5 


-29±12 






c2 


0.06±0.G4 


4.11±0.20 


68±2 


37±7 


-31±20 






cl 


0.04±0.03 


3.69±0.20 


70±2 


38±8 


-26±20 




1999/10/06 


Cl 


0.19±0.05 


5.16±0.15 


64±1 


8.9±1.6 


-27±8 






c2 


0.07±0.03 


4.63±0.15 


67±1 


20.0±4.5 


-37±10 






cl 


0.11±0.04 


3.90±0.15 


69±1 


10.9±3.4 


-31±10 




1999/12/05 


Cl 


0.12±0.05 


5.46±0.15 


64±2 


14.5±3.6 


-28±10 






c2 


0.06±0.03 


4.98±0.20 


66±2 


27±5 


-43±15 






cl 


0.10±0.04 


4.11±0.20 


68±2 


18.0±4.0 


-50±15 




2000/04/05 


c2 


0.05±0.03 


5.37±0.20 


64±2 


22±5 


-23±20 






cl 


0.09±0.04 


4.52±0.20 


64±2 


15.6±4.5 


-42±20 




2000/07/17 


c2 


0.05±0.03 


5.90±0.20 


64±2 


22.0±4.5 


-23±20 






cl 


0.09±0.04 


4.69±0.20 


64±2 


13.3±4.0 


-42±20 



Note. — 



The table appears in full form electronically. 
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Table 6. Total and Polarized Flux Spectral Indices 



Source 

(1) 


Group 

(2) 


^opt/lmm 


^opi/ 1mm 




(^>) 






aP 

^mm 

((i) 




3C 111 


LVP 


0, 


.55+0, 


.02 





.60+0, 


.04 


0, 


.52+0, 


.05 


0. 


,15+0. 


,20 


3C 120 


LVP 


0, 


.56+0, 


.02 


1 


.03+0, 


,05 


0, 


.47+0, 


.02 


-0, 


,28+0, 


,34 


3C 273 


LVP 


0, 


.66+0, 


.01 


1 


.04+0, 


,07 


0, 


.64+0, 


.16 


0, 


,51+0. 


,27 


3C 66A 


IVP 


0, 


.60+0, 


.01 





.50+0, 


,02 


0, 


.16+0, 


.09 


-0, 


,11+0, 


,21 


3C 279 


IVP 


1, 


.19+0, 


.01 


1 


.01+0, 


.04 


0, 


.08+0, 


.06 


-0, 


,13+0, 


.03 


3C 345 


IVP 


1, 


.06+0, 


.01 





.95+0, 


.01 


0, 


.53+0, 


.06 


0. 


,06+0. 


,36 


1803+784 


IVP 














0, 


.28+0, 


.17 


0. 


,19+0. 


,15 


1823+568 


IVP 














0, 


.27+0, 


.17 


0. 


,12+0. 


,23 


BL Lac 


IVP 


0, 


.50+0, 


.01 





.54+0, 


.02 


0, 


.15+0, 


.04 


0. 


,08+0. 


,09 


0420-014 


HVP 


1, 


.00+0, 


.01 


0, 


.86+0, 


,01 


0, 


.25+0, 


.03 


0. 


,11+0, 


,05 


0528+134 


HVP 


0, 


.96+0, 


.02 





.98+0, 


,05 


0, 


.48+0, 


.03 


-0, 


.32+0, 


,26 


OJ 287 


HVP 


0, 


.90+0, 


.01 


0, 


.80+0. 


.01 


0, 


.24+0, 


.06 


-0, 


.25+0, 


.12 


1510-089 


HVP 


0, 


.90+0, 


.01 


1, 


.13+0. 


.04 


0, 


.68+0, 


.07 


-0. 


.44+0. 


.17 


CTA 102 


HVP 


1, 


.00+0, 


.02 


0, 


.97+0, 


,06 


0, 


.37+0, 


.20 


0. 


,11+0. 


,12 


3C 454.3 


HVP 


0, 


.97+0, 


.01 


1 


.05+0, 


,02 


0, 


.40+0, 


.13 


-0, 


,11+0, 


,32 



Note. — uncertainties for aopt/imm are determined from uncertainties of flux measurements; uncertainties 
for amm are based on consistency between different wavelengtlis 



Table 7. Rotation Measures 



Source 


z 


Group 


AxSmm 


Ax7mm 






^^^8GHz 


a 








deg. 


deg. 


lO^rad m-2 


lO^rad m-2 


lO^rad m-2 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


3C 66A 


0.444 


IVP 


-3 


-16 


12+6 








3C 279 


0.538 


IVP 


-3 


-16 


14+2 


-5.7 (2) 


3.1 (1) 


0.92+0.04 


3C 345 


0.595 


IVP 


-2 


-9 


11+7 




-0.33* (3) 


1.6 


1803+784 


0.68 


IVP 


-3 


-16 


12+9 




0.56 (1) 


1.8 


1823+568 


0.664 


IVP 


1 


7 


-10+1 


1.1 (2) 


0.36 (1) 


1.81+0.02 


BL Lac 


0.069 


IVP 


-3 


-16 


6.7+5.5 


7.0 (2) 


0.42 (1) 


1.5+1.2 


0420-014 


0.915 


HVP 


4 


19 


-24+3 








0528+134 


2.06 


HVP 


-16 


-87 


280+30 




1.5 (1) 


3.1 


OJ 287 


0.306 


HVP 


6 


30 


-19+11 








1510-089 


0.361 


HVP 


-7 


-38 


26+10 








CTA 102 


1.037 


HVP 


18 


97 


-140+40 




-4.2* (4) 


1.6 


3C 454.3 


0.859 


HVP 


-4 


-19 


24+13 




9.0 (1) 


1.8 



Note. — Columns are as follows: (1) - name of the source; (2) - redshift; (3) - polarization variability 
group; (4) - RM correction at 3 mm; (5)- RM correction at 7 mm, (6) - intrinsic rotation measure obtained 
in the 43 GHz core; (7) - intrinsic rotation measure obtained in the 15 GHz core; (8) - intrinsic rotation 
measure obtained in the 8 GHz core; * - intrinsic rotation measure obtained in the core at 5 GHz instead of 



GHz; nu mbers in parentheses de note referenc e s for R M° measured at lo w frequencies: 1 - IZavala fc Tavlor 



(l2003h . 2 

RM (X d 



Gabuzda et al. I ([20061), 3 - iTavlor I (|l998[ ). 4 - iTavlor I (|2000[ ): (9) - the exponent in the relation 
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Table 8. Correlation Coefficients between Polarization in the Core at 7 mm 
and Overall Polarization at Optical, 1 mm, and 3 mm Wavelengths 



Source 


.fr,7 


/fr 


J 37 




rT7 


'37 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


3C 66A 


0.81 (4) 


0.22 ( 8) 


0.81 (3) 


0.97 (4) 


0.21 (8) 


0.81 (3) 


3C 111 


••• (1) 


... (0) 


... (2) 


... (0) 


... (0) 


... (0) 










n 79 


v.oo y 


{I) 


3C 120 


... (1) 


... (2) 


... (1) 


... (1) 


... (2) 


... (1) 


0528+134 


... (1) 


-0.17 ( 8) 


-0.38 ( 3) 


... (1) 


-0.49 ( 8) 


0.14 (3) 


OJ 287 


0.33 (6) 


0.78 ( 7) 


-0.11 (5) 


0.30 (6) 


0.68 ( 7) 


0.25 (5) 


3C 273 


... (1) 


... (2) 


... (2) 


... (1) 


... (2) 


... (2) 


3C 279 


0.31 (5) 


0.34 ( 9) 


0.85 (6) 


0.21 (5) 


0.39 ( 9) 


0.90 (6) 


1510-089 


... (1) 


0.79 ( 5) 


0.97 (4) 


0.08 (3) 


0.65 ( 5) 


0.92 (4) 


3C 345 


-0.19 (4) 


0.30 ( 6) 


0.59 (5) 


-0.03 (5) 


0.79 ( 6) 


0.36 (5) 


1803+784 


... (0) 


-0.09 ( 5) 


-0.09 (5) 


... (0) 


-0.70 ( 5) 


0.61 (5) 


1823+568 


... (0) 


0.29 ( 5) 


0.76 (4) 


-0.89 (3) 


0.30 ( 6) 


0.94 (4) 


BL Lac 


-0.97 (3) 


-0.16 (12) 


0.98 (3) 


0.18 (4) 


0.23 (12) 


0.81 (4) 


CTA 102 


-0.68 (3) 


0.34 (11) 


0.90 (4) 


-0.78 (3) 


0.18 (12) 


0.72 (4) 


3C 454.3 


... (2) 


-0.17 (11) 


-0.14 (3) 


-0.41 (4) 


-0.14 (11) 


-0.59 (3) 



Note. — In this and subsequent tables, coefficients of correlation indicated by bold font are significant at 
a significance level e=0.1; integers in parentheses indicate number of observations. 



Table 9. Correlation Coefficients between Polarization in the Inner Jet at 
7 mm and Overall Polarization at 1 mm and 3 mm 



Source 


F'l7 


F'37 




^ 37 


(1) 


(2) 


(3) 


(4) 


(5) 


3C 111 


0.25 ( 7) 


... (2) 


0.64 ( 7) 


... (2) 


3C 120 


0.87 ( 5) 


... (1) 


0.78 ( 6) 


... (1) 


3C 273 


0.79 ( 8) 


0.98 (5) 


0.41 ( 9) 


0.99 (5) 


CTA 102 


0.70 (12) 


0.51 (4) 


0.53 (12) 


0.80 (4) 


3C 454.3 


0.36 (13) 


0.22 (3) 


0.61 (13) 


0.12 (3) 
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Table 10. Correlation Coefficients between Degree of Polarization and Total 

Flux Density 



Source 


^opt 




^3 mm 




(1) 


(2) 


(3) 


(4) 


(5) 




n ns (P.\ 
u.uo yo ) 




U.O / l^lU ) 




"^P 111 




0^1 {^ ^^'\ 

— U.Ol V J 


U.Uu \ Oj 


... ( n\ 


0/190 01/1 




c\f\ ^^ ^\ 

— U.UD ) 


( 1 ^ 
{ I) 




3C 120 


■■■ (2) 


-0.48 ( 8) 


... (2) 


... (2) 


0528+134 


... (1) 


-0.33 (11) 


-0.89 ( 3) 


0.02 (17) 


OJ 287 


0.10 (8) 


0.02 (17) 


-0.36 (13) 


-0.28 (17) 


3C 273 


-0.01 (6) 


0.20 (20) 


0.32 (20) 


-0.06 ( 5) 


3C 279 


-0.13 (7) 


0.23 (22) 


-0.13 (18) 


0.13 (17) 


1510-089 


... (2) 


-0.62 ( 5) 


0.74 ( 9) 


-0.18 (17) 


op OA C 


0.98 (5) 


-0.41 (16) 


0.34 (22) 


-0.23 (17) 


1803+784 


... (0) 


0.11 ( 6) 


-0.35 (20) 


0.48 (17) 


1823+568 


... (0) 


0.00 ( 8) 


0.43 (20) 


0.38 (17) 


BL Lac 


-0.57 (5) 


-0.51 (19) 


-0.07 ( 6) 


-0.12 (17) 


CTA 102 


-0.43 (3) 


0.25 (13) 


0.36 ( 5) 


0.20 (17) 


3C 454.3 


0.44 (3) 


0.14 (14) 


-0.19 ( 5) 


0.67 (17) 


Table 11. 


Correlation Coefficients between Polarized and Total Flux Density 


Source 


/opt 


^Imm 




/7mm 


(1) 


(2) 


(3) 


(4) 


(5) 




0.61 (6) 


0.12 (11) 


0.60 (10) 


0.82 (17) 




... (1) 


0.24 (10) 


0.15 ( 3) 


... (0) 


O/i on m A 
U4ZU — Ui4 


0.92 (5) 


0.88 (11) 


... (1) 


0.84 (17) 


3C 120 


... (2) 


0.05 ( 8) 


... (2) 


... (2) 


0528+134 


... (1) 


0.32 (11) 


-0.81 ( 3) 


0.55 (17) 


OJ 287 


0.72 (8) 


0.44 (17) 


0.16 (13) 


0.60 (17) 


3C 273 


0.21 (6) 


0.76 (20) 


0.57 (20) 


0.49 ( 5) 


3C 279 


0.78 (7) 


0.60 (22) 


0.61 (18) 


0.69 (17) 


1510-089 


... (2) 


0.03 ( 5) 


0.81 ( 9) 


0.06 (17) 


3C 345 


0.99 (5) 


0.21 (16) 


0.74 (22) 


0.39 (17) 


1803+784 


... (0) 


0.26 ( 6) 


-0.20 (20) 


0.66 (17) 


1823+568 


... (0) 


0.77 ( 8) 


0.72 (20) 


0.79 (17) 


BL Lac 


-0.14 (5) 


0.30 (19) 


0.67 ( 6) 


0.66 (17) 


CTA 102 


-0.30 (3) 


0.60 (13) 


0.52 ( 5) 


0.76 (17) 


3C 454.3 


0.63 (3) 


0.80 (20) 


0.08 ( 5) 


0.91 (17) 
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Table 12. Average Values of Fractional Polarization 



Group 




Opt 


1 mm 


3 mm 


7 mm 


LVF 


/ \ ro/i 


1.4±1.2 (1.0) 


10.2±2.1 (7.2) 


2.8±0.9 (2.0) 


1.4±0.3 (1.0) 




jyinax 
cell 


3104 


47 


625 


2500 


HVP 


(Wmax) [%] 


11.4±7.5 (8.1) 


8.6±1.8 (6.1) 


3.1±1.3 (2.2) 


6.0±1.9 (4.2) 




Armax 
-'^cell 


47 


66 


510 


136 


IVP 


(Wmax) [%] 


29.0±11.3 (20) 


16.8±9.4 (11) 


9.1±2.4 (6.4) 


7.5±2.1 (5.3) 




Armax 
■'^cell 


8 


18 


59 


87 


LVP 




0.4±0.3 (0.3) 


1.3±0.3 (0.9) 


0.9±0.6 (0.6) 


O.liO.l (0.1) 






38025 


2900 


6050 


490000 


HVP 




1.7±1.2 (1.2) 


1.6±0.8 (1.1) 


0.6±0.3 (0.4) 


0.5±0.2 (0.3) 






2105 


1914 


13611 


19600 


IVP 


(m„,i„)[%] 


10.6±8.2 (7.5) 


4.0±3.1 (2.8) 


3.2±2.1 (2.3) 


2.0±1.0 (1.4) 




^SiT 


54 


306 


480 


1225 



Note. Values in parentheses indicate standard deviations expected from the cellular model with number 
of cells derived from the fractional polarization; see text 



Table 13. Derived Shock Strength 



Source 




»?opt 




^Imm 




^3min 




^7mm 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


3C 66A 


29.7 


1.94 


36.4 


3.20 


9.3 


1.21 


7.7 


1.17 


116.0 


3C 111 


3.1 


1.05 


12.2 


1.22 


2.3 


1.04 


0.5 


1.01 


108.3 


0420-014 


26.2 


1.81 


7.5 


1.18 


0.5 


1.01 


3.3 


1.07 


59.8 


3C 120 


0.5 


1.01 


11.1 


1.30 


2.1 


1.05 


1.1 


1.02 


124.5 


0528+134 


7.0 


1.15 


12.1 


1.32 


1.4 


1.03 


5.4 


1.12 


57.0 


OJ 287 


15.6 


1.24 


6.9 


1.11 


4.0 


1.06 


5.0 


1.08 


85.3 


3C 273 


0.5 


1.01 


7.2 


1.11 


4.0 


1.06 


1.6 


1.02 


96.8 


3C 279 


39.2 


3.49 


10.9 


1.30 


13.7 


1.40 


6.9 


1.17 


59.2 


1510-089 


4.1 


1.06 


9.7 


1.15 


4.8 


1.07 


9.7 


1.15 


83.6 


3C 345 


38.3 


1.78 


11.4 


1.18 


8.1 


1.13 


9.6 


1.15 


82.7 


1803+784 






11.3 


1.24 


5.9 


1.11 


5.7 


1.11 


67.1 


1823+568 


29.8 




20.3 




9.9 


1.91 


10.6 


2.08 


37.0 


BL Lac 


7.8 


1.12 


10.4 


1.18 


7.8 


1.13 


4.9 


1.08 


86.3 


CTA 102 


9.4 


1.14 


7.1 


1.12 


3.4 


1.06 


6.1 


1.10 


75.9 


3C 454.3 


6.2 


1.28 


8.5 


1.51 


1.8 


1.07 


6.6 


1.35 


39.0 
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Fig. 1. — Degree of polarization at optical, 1 mm, and 3 mm wavelengths from the whole 
source vs. degree of polarization measured in the VLBI core at 7 mm for the most nearly 
simultaneous pair of observations for each source. The linear coefficient of correlation, r, 
is given in each panel. Symbols denote the quasars (open circles), BL Lac objects (filled 
circles), and radio galaxies (triangles). 
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Fig. 2. — Left panel: Dependence between polarization variability index at optical, 1 mm, 
and 3 mm wavelengths and polarization variability index in the VLBI core at 7 mm. Right 
panel: Dependence between polarization position angle variability index at optical, 1 mm, 
and 3 mm wavelengths and polarization position angle variability index in the VLBI core at 
7 mm. The symbols are the same as in Fig. [H 
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Fig. 3. — Connection between polarization and position angle variability indices in the VLBI 
core at 7 mm. Symbols denote the quasars (open circles), BL Lac objects (filled circles), and 
radio galaxies (triangles) 



-so- 




lo 20 30 40 10 20 30 40 10 20 30 40 

Polarization (%) 



Fig. 4. — Distributions of degree of linear polarization in the IVP (left), HVP {middle), and 
LVP {right) variability groups at optical, 1 mm, and 3 mm wavelengths and in the VLBI 
core at 7 mm. 
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Fig. 5. — Total (filled circles) and polarized (filled triangles) fiux density spectra and polar- 
ization position angle (crosses) measurements from optical to 7 mm wavelengths. 
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Fig. 6. — Dependence between polarized and total flux spectral indices calculated between 
the optical and 1 mm wavelengths (open circles - HVP sources, filled circles - IVP sources, 
triangles - LVP sources). 
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Fig. 7. — Dependence of polarization position angle on square of wavelength in the IVP 
and HVP sources. Sohd line represents approximation of the dependence by a Faraday 
rotation law. 




Fig. 8. — Distributions of offsets between EVPAs at 3 mm and in the 7mm core before the 
RM correction {left panel) and after the RM correction [right panel). The distributions do 
not include the observations used to calculate the RM. 



-55- 




Fig. 9. — Dependence of polarization position angle on square of wavelength in 3C 273. Solid 
line represents approximation of the dependence by a Faraday rotation law. 
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Fig. 10. — Fractional polarization vs. total flux density at 1 mm from the whole source and 
in the core at 7 mm. Filled circles correspo nd to our observations and triangles repesent the 
1 . 1 mm from iNartallo et al.l (119981 ) . 
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Fig. 11. — Distributions of offsets between EVPAs for simultaneous observations at different 
wavelengths in the IVP and HVP sources. 
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Fig. 12. — Distributions of offsets between EVPAs for simultaneous observations at different 
wavelengths in the LVP. The polarization position angle at 7 mm corresponds to the EVPA 
in the inner jet. 
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Fig. 13. — Distributions of offsets between polarization position angle and jet direction. 
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Fig. 14. — Dependence of degree of polarization on wavelength in the IVP (left panel) 
and HVP (right panel) sources. Solid lines represent the expected Faraday depolarization 
according to eq. H] and RM derived in §5.1. In general, the fit is poor, so that in most of 
our sources Faraday depolarization is unable to explain the frequency dependence of the 
polarization. 
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Fig. 15. — Lorentz factor vs. polarization variability index at different wavelengths. Symbols 
have the following meaning: triangles - LVP sources, filled circles - IVP sources, open circles 
- HVP sources. 
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Fig. 16. — Sketch of the main emission regions in the compact jet compatible with the results 
of our investigation. Note the logarithmic scale of distance from the black hole (denoted by 
the dot at the center of the flared accretion disk). The magnetic fleld is a tight helix up 
to the "true " core at 1 mm, beyond which it is chaotic. The system is symmetric with 
respect to the accretion disk (only one side is shown). The spectrum of the ambient jet 
rises gradually up to the peak frequency, at which the optical depth at the outer edge of the 
acceleration zone is of order unity, and cuts off at infrared frequencies. The spectrum of a 
shock moving through the core is more highly inverted below the peak and its high-frequency 
cutoff extends beyond the optical region. Although only one is shown, multiple shocks are 
often present at any given time. 
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Fig. 17. — Left panel: Distributions of offsets between polarization position angle at 1 mm 
and inner jet at minimum flux observed at 1 mm in each source. Middle panel: Distributions 
of offsets between polarization position angle at 7 mm and local jet direction for superlu- 
minal knots listed in Table [51 Right panel: Distributions of the degree of polarization in 
superluminal knots (one entry for each knot; see text). 



